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ABSTRACT 
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on the Catani-Seymour dipole factorisation, as first suggested by [TP] , is 
discussed. First results obtained with the new algorithm are compared with 
experimental data. 



st ef f enOtheory . phy . tu-dresden . de 
2 f rank . kraussOdurham .ac.uk 



Contents 



Introduction 2 

1.1 Parton showers based on subtraction methods [3 

1.2 Short review of the Catani-Seymour subtraction method 

Construction of the algorithm 

2.1 Colour factors and spectators 

2.2 Ordering parameter 

2.3 Scales to be chosen 

2.4 General considerations on massive particles 

Kinematics of the individual splittings 

3.1 Final-state emitter and final-state spectator 

3.1.1 Massive case 

3.1.2 Massless case 

3.2 Final-state emitter and initial-state spectator 

3.2.1 Massive case 

3.2.2 Massless case 

3.3 Initial-state emitter and final-state spectator 

3.3.1 Massive case 

3.3.2 Massless case 

3.4 Initial-state emitter and initial-state spectator 

3.5 SUSY QCD splitting functions 

Comparing the hardest emission with matrix elements 

4.1 Three-jet production at lepton-colliders 

4.2 Real corrections to leading order DIS 

4.2.1 The gluon emission process 

4.2.2 The initial-state gluon channel 

4.3 Associated production of a weak gauge boson and a light parton 

4.3.1 The gluon emission channel 

4.3.2 The initial-state gluon case 

Applications 

5.1 Jet production at e + e~ colliders 

5.1.1 Comparison with LEP1 data 

5.1.2 Jet rates in heavy-quark production 

5.2 Particle production in hadron collisions 

5.2.1 Inclusive gauge boson production 

5.2.2 Inclusive jet production 

Conclusions and outlook 



1 



1 Introduction 



In the past decades, parton shower Monte Carlo programs, such as PYTHIA [3lH] or HERWIG 
[5.6J have been indispensable tools for planning and analysing particle physics experiments at 
different colliders. It can be anticipated that they will play a similarly prominent role in the 
LHC era. 

There are a number of reasons for the success of these workhorses. One of the most important 
ones rests in their ability to bridge the gap between few-parton final states, as described by 
fixed-order perturbative calculations, and the real world, where a multitude of hadrons etc. fills 
the detectors of the experiments. The transformation of the partons of perturbation theory 
into the visible hadrons, hadronisation, is a direct consequence of the confinement property of 
QCD. Currently, this phenomenon can be described in terms of phenomenological models only, 
which depend on various phenomenological parameters tuned to data. These parameters and 
hence the validity of the models in turn depends on the properties (such as the flow of energy 
and other QCD quantum numbers) of the parton ensemble; therefore it is important that these 
properties are kept under control. It is the merit of parton showers that they provide a well- 
understood, theoretically sound and universal framework of translating the few-parton states of 
fixed-order perturbation theory, calculated at some high scale, with multi-parton states at much 
lower scales, of the order of a few Aqcd, where hadronisation sets in. In so doing the parton 
showers help guarantee the validity of the tuned parameters of the hadronisation models. 

To achieve this translation of few-parton to multi-parton states, the parton shower programs 
rely on correctly describing QCD particle production in the dominant soft and collinear regions 
of phase space, giving rise to the bulk of radiation. It is in this region, where the complicated 
radiation pattern of multiple particle emission factorises into nearly independent - up to or- 
dering in terms of a suitably chosen parameter - individual emissions of single partons. This 
approximation, namely expanding around the soft or collinear limit, ultimately leads to the 
resummation of the corresponding leading logarithms, which are then typically encapsulated in 
exponential form in the Sudakov form factors. Their probabilistic interpretation in fact is the 
central feature allowing for a straightforward implementation in an event generator, producing 
unweighted events. Due to the resummation of leading logarithms it should thus not be too 
surprising that the parton shower programs more than often produce answers for QCD-related 
questions, which approximate exact results very well. 

However, the quality of the answers provided by the parton shower approximation alone 
relies on whether the question is related to the soft and/or collinear region in the phase space 
of particle production. If this is not the case, for example because of the relevance of hard 
emissions or of non-trivial correlations of particles, the quality of parton shower results tends 
to deteriorate. In such cases, evidently a full quantum mechanical treatment as provided by 
fixed-order calculations becomes mandatory. Therefore, the problem of systematically including 
higher-order effects into parton shower programs has been in the center of research in the 
past few years. In principle, there have been two major avenues of investigation. One dealt 
with the question of how to include the correct QCD next-to leading order correction to total 
cross sections [Zj-[9], and has led to an implementation ready for use by the experiments in 
form of the MC@NLO code [TU]. The other considered the inclusion of tree-level multi-leg 
matrix elements into the simulation [H]-[I3], and has lead to two types of algorithms being 
implemented. One, based on pnfT2] . is the cornerstone of the event generator SHERPA [16j-[18j 
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and an alternative formulation of the same algorithm, proposed in [13], has been implemented 
in ARIADNE [T9.20J. The other merging algorithm, based on [T^lTTo] . has been incorporated in 
ALPGEN [2TJ, MAdGraph/MadEvent [22]- [23, and HELAC fI5TIU\ . Although it is not entirely 
clear how these two approaches relate in detail, some first comparisons [27] show an interesting 
and assuring degree of agreement. 

As one of the most recent outcomes of this line of research, it became apparent that in order 
to systematically improve the event generators by including higher-order corrections, also the 
parton shower algorithms themselves must be ameliorated. Some developments in this direction 
include an improved treatment of angular ordering and massive partons in HERWIG++ [28] or 
the introduction of a new fcj_-ordered shower in PYTHIA [29J. More recently, and motivated by 
the wish to include loop-level calculations in a more straightforward and systematic manner, the 
application of subtraction terms, prevalent in QCD next-to leading order calculations, has been 
proposed. This paper reports on the construction of a parton shower based on such subtraction 
terms. It uses the Catani-Seymour dipole formalism [3T)f3"T] and the corresponding subtractions 
as a starting point 0. This formulation of a parton shower has been proposed for the first time 
in [Tf2] . A similar ansatz relies on antenna subtraction terms [33f34"] and has been presented 
recently in [35] . 

The paper is organised as follows: After briefly introducing the idea of parton shower algo- 
rithms based on subtraction terms in Sec. 11.11 and a short review of the subtraction formalism 
of Catani-Seymour in Sec. 11.21 Sec. [2] states the basic construction principles of the proposed 
shower description. In Sec. [3] the actual parton shower built on Catani-Seymour subtraction 
terms is constructed. The most general massive and the massless case for all the possible QCD 
splitting types are discussed in detail, and the modifications needed to include splittings of 
super symmetric particles are discussed. The analytic expressions for the first shower emission 
from various core processes are compared with the corresponding exact tree-level matrix ele- 
ment calculations in Sec. HI In Sec. [5] predictions obtained with the developed shower formalism 
are confronted with experimental data and other calculations. The focus hereby is on hadron 
production in e + e~ collisions, and Drell-Yan and QCD jet production at the Fermilab Tevatron. 
Sec. E] is devoted to the summary, conclusions, and an outlook on further developments. 

1.1 Parton showers based on subtraction methods 

Since its formulation almost a decade ago, the Catani-Seymour dipole formalism [3T)1I3"T] has 
been widely used in the calculation of next-to leading order (NLO) corrections in QCD, see for 
instance [36]-[i3]. 

Such calculations typically face the problem of infrared divergences in both the real and the 
virtual parts of the NLO correction. In principle, such divergences are not really a problem, since 
for physically meaningful observables, the Kinoshita-Lee-Nauenberg theorem |44"f4"5j guarantees 
their mutual cancellation. To technically perform this cancellation, however, the divergences 
need to be regularised, which is usually performed by dimensional regularisation, i.e. continuing 
the calculation to d dimensions. There, the infrared divergences manifest themselves as poles 
in 1/(4 — d) or 1/(4 — d) 2 . To deal with the poles and achieve the cancellation, subtraction 
methods may be used. In general, they rely on the fact that the infrared divergences in the 

1 This approach has also been employed in a parallel project, 32 . 
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real correction part follow an universal pattern. This allows to construct simplified terms in a 
process-independent way that encapsulate all infrared divergences occurring in the full matrix 
element. Then, subtracting these terms from the real-correction matrix elements will yield 
an infrared-finite result, such that this subtracted matrix element can be safely integrated 
numerically in four dimensions. In addition, the subtraction terms are chosen such that they 
can be analytically integrated in d dimensions over the phase space of the additional soft or 
collinear particle causing the divergences. This yields the poles in 1/(4 — d) or 1/(4 — d) 2 , which 
are then added to the virtual part of the correction, and thus cancel the poles there. 

The catch with the subtraction methods is that the subtraction terms can be constructed 
locally from the (colour-ordered) Born matrix element. In the Catani-Seymour method, for 
instance, pairs of particles are interpreted as emitting particle and spectator and are subjected 
to a splitting kernel creating a third particle. In this splitting process, one of the particles 
actually splits, while the recoil is compensated for by the spectator, which may be interpreted 
as its colour partner. At the same time, the phase space factorises exactly into a phase space 
over the original particles, already present at the Born level, and into a phase space of the 
additional particle emerging in the splitting. This exact factorisation corresponds to an exact 
mapping of the two original momenta (emitter and spectator) onto three four-momenta. At 
each point of the procedure all particles remain on their respective mass shell. 

This is why constructing parton showers based on such methods currently is being pursued 
by different groups. It is clear that these showers, in full conformance with original formula- 
tions employing the splitting of individual, single partons, are based on the universal soft and 
collinear dominance of QCD radiation. Similar to the original shower algorithms, the emerg- 
ing large logarithms occurring with each individual parton emission can be resummed in a 
straightforward way through a Markovian process. This, in principle, renders both formula- 
tions formally equivalent. On the other hand, however, showers based on subtraction terms 
have the practical advantage that the conservation of four-momentum is built in with particles 
that remain on their mass shell at any given point H. It can be anticipated that these features 
ultimately will allow for a more transparent merging with multi-leg matrix elements and a 
drastically alleviated matching with full NLO calculations. 

1.2 Short review of the Catani-Seymour subtraction method 

The Catani-Seymour subtraction method has been introduced in [3U] for massless partons and 
it has been extended to massive partons in [3T]. To fix the notation for the rest of the paper, 
it will be briefly reviewed here. 

The essence of this method is embedded in the dipole factorisation formula 

\M m+1 r = EE + E E n + E E v t + EE ^ ai ' b + ■■■■ (i) 

i,j kj^i,j i,j a a,i k^i a,i b^=a 

The individual dipole contributions D provide the correct approximation of the (m + l)-parton 
matrix element squared in the different singular regions of phase spaced. In each term i, j and 

2 It is interesting to note that the latest refinements of the parton showers in HERWIG and PYTHIA also put 
more emphasis on the notion of a colour-connected partner compensating recoils etc. [28|29j . 

3 Note that squared matrix elements shall always be understood as properly normalised with respect to the 
colour degrees of freedom of incoming particles. 
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k denote final-state partons and a and b stand for initial-state partons. The first sum always 
runs over the two particles to be combined, whereas the second sum takes care of the spectators. 
Accordingly, the four terms correspond to the splitting of a final-state parton accompanied by 
a final-state or initial-state spectator and emissions off incoming particles in the presence of 
a final-state or an initial-state spectator, respectively. Finally, the dots in the equation above 
denote some potential finite terms which do not exhibit any divergence. 

For the case of final-state emitters with a final-state spectator, for instance, the individual 
dipole contributions read 



V 



ij,k 



»<1, 



k. 



m + 1| 



T, ■ T, 



~^ij,k 1 1 



k. 



,m + l) m , (2) 



when all the involved partons are assumed to be massless. The occurring m-parton states are 
constructed from the original (m + l)-particle matrix element by replacing the partons i and 
j with the new parton ij, the emitter, and the original parton k with k, the spectator. In the 
massless case, their momenta are given by 



ij 



Pi 



Vij,k 



1 



Vij,k 



Pk 



and 



Pk 



1 



1 



Vij,k 



Pk 



where the dimensionless, Lorentz-invariant quantity is given by 

PiPj 



yij,k 



PiPj + PiPk + PjPk 



(3) 



(4) 



It is simple to show exact four-momentum conservation, i.e. pfj + p^ = Pi + Pj + Pk> with all 
particles on their mass shell. In the matrix element on the right hand side of Eq. ([2]), the Ty, 
Tfe are the colour charges of the emitter and spectator, respectively, and the Vy ^ are matrices 
in the emitter's spin and colour space, responsible for its branching. The operators "Vij t k also 
depend on the dimensionless, Lorentz-invariant quantities 



PiPk 



PiPk 



PiPk + PjPk PijPk 



and 



PjPk 



= PjPk 

PiPk + PjPk PijPk 



(5) 



For instance, for the case of a quark splitting in the final state with a final-state spectator, i.e. 
Qij — » <?i + 9j, where s and s' denote the spins of ij and i, respectively, and where the subscripts 
label the momenta, 



(s|V 



H9j 



,k[Zi 



■ yij,k)\s') = 87r/i e a s C F 



Hi 



Vij,k) 



(6) 



Here, e = 4 — d, with d the number of dimensions. Similar expressions emerge for the other 
QCD splittings or when masses are included. However, as a general property, the matrices 
do not become singular, if any of the scalar products PiPj, PiPk or pjp k vanishes, and therefore 
the only soft or collinear divergences in the dipole terms Py,fe are related to p^pj — > 0. 

The collinear limit of the two final-state partons i and j originating from a splitting ij — > i+j 
is defined through their relative transverse momentum k± — > 0. This limit can be investigated 
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by decomposing the momenta as 



£*2 it fJ> 

rf = ^ + (7) 

= d-.y + ^-*i, (8) 

where the lightlike p^ defines the collinear direction and is an auxiliary lightlike vector that 
specifies the spacelike transverse momentum k±, with k\ = — k^_, through pk± = nk± = 0. 
Then, in the collinear limit, the scalar product PiPj reads 



and the dipole variables are given by 



Vij,k 



kl 



Zi 



2z(l - z)pp k ' ~ l ' ~ 3 
Pk Pk and Pij^P^- (10) 

It can then be shown that in this limit the matrices become proportional to the Altarelli- 
Parisi splitting kernels, 

V^fc -> 87T/i 2e a; s P(ij),i(z, k ± ; e) . (11) 

In this limit the only remaining dependence of the dipole contributions 'Dij k on t ne spectator 
fc resides in its colour factor and it can be shown that Eq. (j2J) reproduces the well-known 
universal collinear behaviour of the (m + l)-parton matrix element, 

m+i(l, ...,i,...,j,...,m+ 1||1, j, ...,m+ l) m +i 

kj^o 4tt/i a s ij m + llPujuiz, k ± ; e)|l, . . . , ij, . . . , m + l) m , (12) 

PiPj 

where again, the kernel P is a <i-dimensional Altarelli-Parisi splitting function. 

In contrast, the limit where pj becomes soft is given by pj = Ag M with A — > and g M some, 
in principle arbitrary, four- vector. In this limit, the dipole variables become 

Uij,k * , Z{ 1 Zj > 1 , 

Pk Pk and Pij (13) 

and Vj^fc tends to 

1 A ->0 1 PiPfc 



z~i(l-y ij)k ) X (pi+Pk)q 



(14) 



Therefore, 



It can thus be shown that the well-known soft limit of the (m + l)-parton matrix element is 
recovered, namely 

m+i(l, j, ■ ■ ■ , m + 1||1, . . ., i, j, ■ . ■ , m + l) m+ i 

A ^o \- 87r/x 2e g s .. T fc -T^(p^) .. 

-< 1 .-.«.-.™+ H H,...,«,...,m + 1>-.(W) 

Taken together, these considerations and similar reasoning for the other dipole contributions 
translate into the dipole formula, Eq. (1T1), to provide a point-wise approximation to the full 
(m + l)-parton matrix element, which exactly recovers all the soft and collinear divergences. 

Before starting the discussion on the construction of a parton shower algorithm from the 
Catani-Seymour dipole formula in Sec. [2] the generalisation of Eq. ([9]) to the massive case and 
the analogous result for the splitting of an initial-state parton shall be briefly repeated. 

First, re-consider the splitting ij —>■ i + j from above. This time, however, both the emitter 
and the splitting products are allowed to be massive, the corresponding mass shell conditions 
read p 2 = m 2 j, p 2 = mf and p 2 = m 2 . The momenta pi and pj can again be written in a 
Sudakov parametrisation according to 

Pi = ZP^ + " 7^ + k l, ( 17 ) 

z 2pn 

, „ -k 2 ± - (1 - z) 2 ml + m 2 

with n 2 = and k± perpendicular to both p and n. Identifying k\ = — h 2 ± the invariant mass 
of partons i and j is now given by 

Accordingly, the collinear singularity is shielded when at least one of the two partons has a 
finite mass. 

Finally, consider the case when final-state parton % becomes collinear to an initial-state 
parton a. This corresponds to the splitting a — > ai + i, with ai the initial-state parton that 
enters the m-parton process. Considering only massless initial states, all the partons involved 
in the splitting are consistently taken to be massless. Decomposing the final-state momentum 
Pi according to 

pt = (l-a;)j£+_-L- + k ± , (20) 

1 — x 2p a n 



the collinear limit is reached for 



™'=m^)' ki ^°- (21) 



2^ 



with the magnitude of the spacelike transverse momentum vector k±, namely k\ = — k 
The definitions Eq. (TITH) and Eq. fl2Tl) constitute the basic relations for identifying the transverse 
momentum vector for the different splitting types in terms of the respective splitting variables 
used to describe the branchings, see Sec. El 
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2 Construction of the algorithm 



To formulate a parton shower algorithm based on the Catani-Seymour dipole formulae, the 
corresponding splitting operators T> that describe the emission of an additional parton from an 
arbitrary m-parton state have to be analyzed and rewritten in a suitable form, before they can 
be used for a showering algorithm. To this end, a number of issues has to be resolved: 

• First of all, only the four- dimensional expressions of the splitting kernels T> will enter 
the parton shower. In addition, the splitting kernels are employed in their spin-averaged 
form. This manipulation is straightforward and a detailed discussion is therefore not 
necessary. The resulting splitting kernels depend on the actual configuration of emitters 
and spectators in the initial- and final state and they will be listed in the corresponding 
parts of Sec. [3j 

• In order to keep the probabilistic notion enabling simulation, to use a Markovian for- 
mulation for the showering process and to facilitate the hadronisation at the end of the 
shower, issues concerning colour correlations have to be solved. While the original Catani- 
Seymour dipole formulae consider all colour correlations, the shower will account only for 
the leading terms in 1/N C . This will be further discussed in Sec. 12.11 

• Also, the phase space factorisation and the corresponding combination procedure is ef- 
fectively inverted to construct the kinematics of the individual splittings. This yields 
splitting kernels for 1 — ► 2 QCD branchings that allow for the inclusion of finite parton 
masses in quite a general way. Each splitting parton thereby is accompanied by a single 
colour-connected spectator parton compensating the recoil of the splitting. The only ex- 
ception here are initial-state splittings in the presence of an initial-state spectator, where 
the recoil is taken by all final-state partons of the event. The introduction of the spectator 
allows to assemble the shower kinematics such that four-momentum conservation can be 
ensured after each individual branching with all external partons on their mass-shell. Ac- 
cordingly, this parton shower algorithm can be stopped at any intermediate stage as well 
as started again for a partially evolved parton ensemble. However, the exact procedure 
for reconstructing the kinematics of each splitting again depends on whether the emitter 
and spectator are in the initial- or final state, respectively. The corresponding formulae 
are listed in Sec. El 

• The actual shower evolution variable specifying and ordering subsequent emissions is 
chosen to be the transverse momentum between the splitting products for branching 
final-state partons and the transverse momentum with respect to the beam for emissions 
from the initial state, collectively denoted by k^. The physics underlying this choice will 
be further detailed in Sec. 12.21 

• Furthermore, choices have to be made concerning the scales entering the QCD running 
coupling constant, a s , and the parton distribution functions when initial-state partons 
are present. This will be discussed in Sec. 12.31 

• Based on these considerations, appropriate Sudakov form factors are constructed that 
determine the probability for a certain branching process not to occur for a given range 
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of the evolution variable, k^. These Sudakov form factors constitute the basis of the actual 
Monte Carlo showering algorithm. Again, their specific form depends on the details of 
emitter and spectator parton and they will thus be given in corresponding parts of Sec. 
El too. 

• This section closes with some general considerations concerning the treatment of parton 
masses, cf. Sec. 12.41 



2.1 Colour factors and spectators 

The starting point for every parton shower evolution is a given set of partons and their momenta 
from a fixed-order matrix element calculation. In the large- N c approximation a colour flow 
can be assigned to each parton configuration. Since in most cases the initial matrix element 
calculation is already summed and averaged over the colours of final and initial partons, the 
assignment typically is performed a posteriori in different ways in different codes. However, as 
a result the partons entering the parton shower after this assignment have a well-defined colour, 
and, due to the large- N c limit, one or two uniquely assigned colour partners @. Motivated by 
considerations on the colour dynamics for soft emissions in the Catani-Seymour formalism, in 
a corresponding shower formulation the spectator parton accompanying a given splitting is 
colour-connected to the emitter parton. For the case of a splitting gluon/gluino then there are 
always two possible colour partners, whereas splitting (anti-)quarks/squarks will have only one 
spectator parton candidate. Following this reasoning, the initial partons will enter the parton 
shower stage in well-defined pairs of potential emitters and spectators. The subsequent parton 
shower will not change this feature. 

To formalise the treatment of colour inside the parton shower presented here, consider the 
colour-operators present in the Catani-Seymour dipole contributions. In the large-iV c limit, 
they are easily calculated for any m-parton state at the price of loosing colour correlations 
beyond 1/ N c . However, in this limit only two cases need to be considered. Independent of the 
actual spectator flavour, the colour algebra for a splitting (anti-)quark/squark yields, 

+ (22) 

whereas a splitting gluon/gluino results in 

+ (23) 

ij * 1 v c 

For convenience, these two results can be combined by introducing Afif ec , the number of possible 
spectators the emitting parton possesses, then 

~^tI~^J^ + ° { nI ) - (24) 



4 Representing the colour flow pictorially by coloured strings of partons, two configurations emerge, namely 
open or closed strings. An open string consists of a colour-triplet state followed by colour octets and ends with a 
colour anti-triplet. Mapping the colour flows, initial-state quarks (colour triplets) correspond to final-state anti- 
quarks (colour anti-triplets), whereas initial-state anti-triplets can be treated as final-state triplets. A closed 
colour string corresponds to a configuration of colour-octet partons only. Accordingly, the end of a closed string 
is colour-connected to its beginning and therefore the whole colour string is invariant under cyclic permutations 
of its individual constituents. 



9 



2.2 Ordering parameter 

Having the individual splitting process under control, i.e. having at hand the corresponding 
splitting kernel with all relevant colour factors and the way the kinematics of the emission 
is constructed, the full showering algorithm with its sequence of splittings can be addressed. 
While the individual splitting kernel properly takes into account the soft and collinear divergent 
regions, in the parton shower itself these regions are cut away and, formally speaking, combined 
with the virtual bits to yield a probabilistic description of the splitting process. The cut on 
the soft and collinear region implies the emergence of corresponding logarithms of the cut 
parameter, which the parton shower aims to resum. Technically, this resummation is achieved 
by arranging the individual emissions in a Markov chain, treating each emission on the same 
footing, and by ordering the emissions with some ordering parameter. This has been detailed in 
textbooks such as [46] . In different parton shower implementations, there are different ordering 
parameters realised, such as the invariant mass of the splitting particle [H]-[IH], the opening 
angle of the pair [281150] . or their relative transverse momentum [T9f2"9] . At the level of doubly 
leading logarithms, these choices are all equivalent, but there are substantial differences on the 
level of next-to leading logarithms, i.e. on the level of single soft logarithms. This is closely 
tied with the treatment of quantum coherence effects [HI]- [HI], which are properly taken into 
account by ordering subsequent emissions through their respective opening angles [50J. In [55J 
it has been shown that another way of properly accounting for coherence effects is evolving in 
a dipole-like picture with subsequent emissions ordered by transverse momenta. 

In the implementation presented here, the parton shower will be ordered by transverse 
momenta, i.e. by the in Eqs. ( IT§]) and fl2Tl) . Apart from the proper treatment of quantum 
coherence effects, this choice has additional benefits: First of all, as will be discussed in the next 
section, cf. Sec. 12.31 by ordering with the ordering parameter also enters as the relevant scale 
in the coupling constant and the parton distribution functions. Second, the definition used here 
allows for a shower formulation on the basis of Lorentz-invariant quantities, see for instance e.g. 
Eqs. OH) and Also, ordering by k^ immediately implies that the parton shower cut-off is 
related to some minimal transverse momentum necessary to resolve partons, which seems quite 
appealing in terms of the physical interpretation of such a resolution criterion. Last but not 
least an ordering by transverse momenta appears to allow for quite a straightforward merging 
of the parton shower with multi-leg tree- level matrix elements in the spirit of [TTp2] . The 
merging method presented there bases on Sudakov suppression weights for matrix elements, 
which are constructed from the transverse momenta of their nodes, and on a vetoed parton 
shower respecting the minimal scale of a fcj_-jet definition. 

In the parton shower evolution each colour- singlet is separately evolved. To this end, all 
emitter-spectator dipoles are iterated over and for each of those configurations a k^ is chosen 
according to the corresponding Sudakov form factor. The dipole with the largest k^ is selected 
to split according to the kinematics detailed below. As long as this largest k^-value is larger 
than the infrared cut-off k_|_ j0 , the shower evolution will continue, and this largest k^ of the 
current evolution step serves as the maximal scale for all dipoles in the colour- singlet in the 
next splitting step. 
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2.3 Scales to be chosen 



When discussing the details of a parton shower implementation, some care has to be taken 
in the choice of various, in principle undetermined, occurring scales. There are a number of 
choices to be made, namely: 

• The evolution variable and the related evolution cut-off: 

As already discussed in the previous section, in this implementation the relative transverse 
momentum of the produced parton w.r.t. its emitter has been chosen as the relevant 
evolution variable. It is given by Eqs. (TK?j) and (I2"T|) . Correspondingly, a cut-off has to be 
set as a tuning parameter, to stay away from phase-space regions where the perturbative 
expansion for the running coupling is divergent. The choice of this cut-off is dictated 
by two aspects. First of all, it seems to be more attractive to try to assign as much 
phase space for particle creation to the, in principle, well-understood perturbative parton 
shower rather than to a phenomenological hadronisation approach such as the Lund string 
fragmentation [56f57] or a cluster model [58]-[60j. This implies that the cut-off should 
be as small as possible. On the other hand, it is clear that perturbative QCD breaks 
down and looses its predictive power at small scales. This is best exemplified by the 
infrared behaviour of the running coupling which exhibits a Landau pole at Aqcd- As 
will be discussed in the next item, since the running coupling in the shower is evaluated 
at a scale related to k^, this feature of QCD prohibits cut-offs in the region of Aqcd- 
Therefore, a suitable choice seems to be a cut-off k^ ]0 of the order of 1 GeV, sufficiently 
separated from the Landau pole. 

• The argument of the running coupling constant, fiR. 

In the previous item it has been already hinted at the choice typically made in parton 
showers, to take the running coupling at scales of the order of k^. The reason for this 
choice is that it incorporates and resums some of the higher-order corrections to the 
splitting. Specifically, in this implementation the choice is to take //^ = [Mr = kj_ if the 
emitter is a final-state particle and fi R s ^ = /j,r = kj_/2 if the emitter is a parton in the 
initial state. 

• The argument of the parton density functions, fip- 

Similar to the case of the running coupling constant, a choice has also been made at which 
scale to take the parton distribution functions, if necessary. In parton showers, there are 
typically two answers, namely to either again take the transverse momentum or to use 
the virtual mass of the initial emitter. Here the choice again is to use fip = k^. 

2.4 General considerations on massive particles 

Taking into account finite quark mass effects in the Standard Model (SM) clearly is of impor- 
tance when producing heavy quarks, bottom or top quarks, in a hard scattering process. In 
addition, many extensions of the SM introduce new strongly-interacting heavy particles, whose 
QCD radiation needs to be modeled to understand the patterns of particle and energy flows in 
their production and eventual decays. Prime examples are scalar quarks and gluinos in super- 
symmetric theories [61] or heavy excitations of the SM quark and gluon fields in models with 
additional space-time dimensions |62j. While at lepton colliders heavy objects only appear in 
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the process' final state, at hadron colliders charm and bottom quarks can also constitute the 
partonic initial state. An example where these are of phenomenological relevance is the asso- 
ciated production of heavy quarks and scalar Higgs particles in supersymmetric models, which 
is a promising channel to gain deeper insight into the mechanism of electroweak symmetry 
breaking, see for instance [63] and references therein. 

In the following section, QCD splitting operators will be derived, that fully take into account 
finite masses of partons in the final state. This includes both emission from heavy particles but 
also the splitting of gluons into heavy quarks such as charm or bottom. Splittings of gluons into 
heavier objects or branchings of heavy states into other heavy objects are beyond the scope 
of this work as they are not well modeled by the soft or quasi-collinear approximation and 
should rather be described with full matrix elements. For all the formulae presented in Sec. [3J 
the massless limit is smoothly obtained when setting the parton masses to zero. This will be 
explicitly examined for some of the important results there. 

Throughout this work, incoming QCD partons will always be treated as massless. The 
leading logarithms that arise for emissions off incoming heavy quarks, logarithms of the type 
(a s \og(Q 2 / rag))™ , with Q 2 the scale of the hard-scattering process and mg the quark mass, 
are summed to all orders in QCD when using heavy-quark parton distribution functions at the 
factorisation scale fip ~ Q and considering the incoming quarks as massless |64,65j. A scheme 
to consistently incorporate explicit masses for incoming heavy quarks, relying on modified 
heavy-quark density functions [66], has recently been presented in [67] . 

3 Kinematics of the individual splittings 

In the following sections, Sees. I3.lti3.4l the actual parton shower built on Catani-Seymour sub- 
traction terms is constructed. To this end, all combinations of initial- and final-state emitter 
and spectator partons are considered in detail, following closely the original publications on 
the subtraction method [30|l3Tj . First, the kinematic variables characterising the individual 
splitting under consideration are discussed. Then the explicit form of the phase space element 
for the three-parton state under consideration is re-expressed through the kinematic variables 
above, and their respective bounds are given. In a next step, the polarisation-averaged split- 
ting kernels for the respective emitter-spectator configuration are listed. This allows to give the 
factorised form of matrix elements with one additional parton in the soft and collinear limits 
of its production and the factorised form of the corresponding differential cross section, which 
includes both matrix element and phase space factorisation. From there, it is quite straightfor- 
ward to deduce the actual Sudakov form factor for the emitter-spectator configuration. Finally, 
the actual kinematics of the splitting is constructed, which may slightly differ from the evolu- 
tion parameters due to mass effects. For each case then also the more familiar massless limit 
is briefly discussed. In Sec. 13.51 the QCD splitting functions for supersymmetric particles are 
presented. 

3.1 Final-state emitter and final-state spectator 

The first case to be investigated is when both the emitter and the spectator parton are in the 
final state, cf. Fig. [T] Accordingly, the splitting {ij,k} — > {i,j,k} has to be studied. When 
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Figure 1: Effective diagram for the splitting of a final-state parton connected to a final-state 
spectator. The blob denotes the m-parton matrix element, and the outgoing lines label the 
final-state partons participating in the splitting. 

considering processes without colour- charged initial-state particles, such as jet production in 
lepton-lepton collisions, this is the only QCD radiation process and thus constitutes the basis of 
a corresponding final-state parton shower. However, the observed factorisation of the differential 
cross section for producing an additional parton also holds in the presence of initial-state 
partons, where only the additional branching channels discussed below then have to be taken 
into account as well. 

3.1.1 Massive case 

In the most general case all partons involved in the splitting can have arbitrary masses, i.e. 
pfj = m?-, pi = p\ = ml, p\ = mf and pj — m|, respectively. In order to avoid on-shell decays, 
which should be described by their respective proper matrix element, only those situations are 
considered, where m?- < mf + m*. 

• Kinematics: 

Exact four-momentum conservation is ensured by the requirement 

Pij +Pk=Pi+Pj+Pk = Q- (25) 

The splitting is characterised by the dimensionless variables y%j^, Z{ and Zj. They are 
given by 

PiPj ~ , ~ PiPk , oft v 

Vij k = , Zi — 1 — Zj — . (26) 

PiPj + PiPk + Pjpk PiPk + Pjpk 

With these definitions the invariant transverse momentum of partons i and j, defined in 
Eq. ( fl9l) . can be written as 

k ± = (Q 2 - -rr?j- ml)y ij)k ^(1 - z { ) - (1 - 5j) 2 m- - z^m] . (27) 

For convenience, the rescaled parton masses 

Tfl 

VQ 
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and the relative velocities between pi + pj and pi (pk), (%,&)> 



- g - fi] - \i\)y\ k - 



[2/il + (1 - $ - /i] - ii 2 k ) (1 - y ij>k )] - A/4 
as well as the velocity between p^ and p k , 



= 1 72~ > ( 31 ) 

are introduced. 
Phase space: 

In the case of a final-state emitter with a final-state spectator, the corresponding three- 
parton phase space d§(pi,pj,p k ; Q) must be analyzed. It exactly factorises into a two- 
parton contribution d§(pij,p k ] Q) and a single-parton phase space factor [dpi(pij,p k )], 

d^(Pi,Pj,Pk, Q) = d$(p ij: p k ; Q) [dpi(p ij: p k )] 6(1 - fi t - ^ - fi k ) , (32) 

where the latter is given by 

[dPi(pij,Pk)] = —, t^r - /== ^ ==^- (1 - Vij,k) dy ijtk dz~ . (33) 

Here and in the following, A denotes the Kallen function, 

X(x, y, z) = x 2 + y 2 + z 2 - 2(xy + xz + yz) . (34) 

The boundaries of the full, unconstrained, phase space read G [0, 2ir], whereas the lower 
and upper limits for z\ and yi^ k are 

y- = — — 2 ' and y+ = 1 - fykix ^ 
i -ni- -fit i - nt - Hj - »l 

respectively. 
Splitting kernels: 

The polarisation-averaged QCD splitting kernels (V^) read 

(v ^ (5 - = * {nbsz " S (1 + z ' + S>} • (37) 

f 1 1 

(v gigj ,k( z i,yij,k)) = 2C mi — - - ~ , + 



1 Zi 4" Zjyij :k Zi + ^fc Ziyij,k 

+ z i (l-z i )-z + z.-2\ 

Vij,k J 

(V QiQ „ fe (^)) = T R — {1-2 [*(!-*)-*+*_]}. (39) 
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Here, Eq. ( 1371) describes the QCD splitting Q — > Qg, of a massive quark Q, the case of a 
splitting anti-quark is formally identical. The corresponding expressions for the splitting 
g — > gg, or g — ► Q<Q are given in Eqs. (j3*81) and (1311 . respectively. Note that in the above 
splitting kernels the free parameter k that occurs in the full NLO subtraction scheme [31] 
has been set to zero to obtain the simplest expressions for the different (V^). 

It should be stressed here that the scalar product piPj present in Eq. ( 1371) can be written 
solely in terms of the splitting variables and the scale k^_: 

Pi p - = -J= + i t—L + — . (40) 

2~ Zi (l-~ Zl ) 2h 2(1-5,) 1 ; 

However, in Eq. ( I37j) the final-state gluon is massless and correspondingly m 2 = such 
that the last term of Eq. ( T4T)1) vanishes in this specific case. 

Matrix element: 

Using the above splitting functions, the full (m + l)-parton matrix element factorises in 
the soft and collinear limit according to 

\M m+1 \ 2 = \M m \ 2 ( T \ _ T T7ip^ 87r «s (Vij,k(Zi>Vij,k)) > ( 41 ) 

ij kjtij \Pi~T~Pj) m ij JV ij 

cf. [30], where the sum covers all the possible emitter-spectator pairs. When combining 
this with the (m + l)-parton phase space a fully factorised expression for the differential 
cross section is obtained, namely 

da m+1 = da m ^ ^ 6% ^ ^ J {y ij>k ) (V ii)fc (£, y ij)k )) , (42) 

a k+n Vat 2n 2lx N a 

where the Jacobian 

emerges from the phase-space factors of Eq. ( l33l) combined with the propagator term of 
Eq. PI}. 

Sudakov form factor: 

A first step toward the construction of the corresponding Sudakov form factor is achieved 
by realising that the ^^-integration in the equation above, Eq. (l4"2l) . can be replaced by 
an integration over the ordering parameter, the transverse momentum, according to 

Uij,k 

Cutting the available phase space through the requirement of a minimal relative transverse 
momentum squared kj_ > k 2 ± > and some upper limit k^ max for the splitting products 



15 



i and j, the Zi integration boundaries become 



i max , kj_ ) 



z + (k] 



k 2 ) 




(45) 
(46) 



with z T taken from Eq. 
be solved for yij,k, 



. Having chosen a valid pair for and z^ this can then easily 



Vij,k 



Q 2 — m 2 — m 2 — ml \Zi(l — z^ 



+ 



1 - Zi)mf 



+ 



ZiTTlj 
1-2 



(47) 



If the calculated fulfils the requirement G [y^,y + ], with y T defined in Eq. 
a valid splitting has been constructed, i.e. a physical branching allowed by phase space. 

The Sudakov form factor corresponding to having no emission from one of the process' 



final-final dipoles between the maximum transverse momentum squared kj 



and the 



infrared cut-off k^_ reads 



A 



FF 



,k I maX ) kj_ ) 



exp 



EE 



V 



\f spec 

ij k^ij yv ij k2 



'dk 2 



2 + 



dzi ^Lkl j ) (y ijk (zi, y ijyk )) 
In 



\ 



(48 



/ 



As already advertised in Sec. 12.31 the scale of the running coupling has thereby been 
chosen equal to the current transverse momentum squared. 

Physical kinematics: 

Having a valid set of splitting variables, the actual physical branching kinematics must 
be constructed in order to fully specify the splitting {ij, k} — > {i,j,k}. In the most 
general case, both the emitter and the spectator parton are massive, prohibiting a simple 
Sudakov parametrisation of Pi and pj in terms of light-like momenta py and pk- Instead 
they must be expressed in light-cone kinematics with massive base momenta. In the 
emitter-spectator centre-of-mass frame the new spectator momentum can be fixed to 



Pk 



2nl + (1 - n 2 - \A - nl)(l - y ijtk ) 



4/4 



+ 



Hi - Hi 



\J H^-j h%-, hVi 

(4)0- - vat) + hI 



v k -- 2 



1 + hI- hI] Q 



Q. 



(49) 



Then the situation is most easily discussed in a frame where Q — pk is at rest and the 
momentum pk points along the z-direction. In this frame, the light-cone momenta of 
Q — Pk and pk can be written as 



Q-Pk = (M,M,0) and p k = (m k e x , m k e x , 0) . 



(50) 
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The ansatz for the light-cone momenta of the new emerging final-state partons reads 

Pi = (rrii,± e y , m it± e~ y , l±) , pj = (m ji± e z , m ji± e~ z , -l±) , (51) 
with m± being the transverse mass of the respective parton, defined according to 



m± = v / m 2 + /|. (52) 

The kinematics is fully determined through energy-momentum conservation and the con- 
straint 

Zi = 1 — Zj = . (53) 

PiPk + PjPk 

Then, 

-* fM 2 + m 2 + m 2 \ 2 /M 2 +m 2 i +m 2 1 -2M 2 z i /coshx\\ 2 

-^-( 2M (imh^jj ' (M) 

and 

M 2 + mf-m 2 coshx / AfsA 

cosh y = — -, smny = — - — cosny , (55) 

2Mm it ± sinni y m^iy 

M 2 -mf + m 2 . coshx / M(l - ^ 

cosh ,2 = — -, smhz = — - — coshz . (56) 

2Mrrij t ± sinni y m j,J- / 

Expressed through ordinary four-vectors the parton momenta in this frame read 

Pi = (m,i t ± cosh y, l± cos 0, l± sin 0, m it ± sinhy) , (57) 
Pj = (m,j t ± cosh z, —l± cos 0, —l± sin 0, m,j t ± sinh z) , (58) 

with the angle not fixed by the splitting and therefore uniformly distributed in the 
transverse plane. The kinematics is completed by rotating and boosting back the momenta 
Pi, pj and pk into the laboratory frame. 

If the spectator is massless, the new final-state momenta can alternatively be given in 
a simple Sudakov parametrisation in the centre-of-mass frame of the emitter and the 
spectator: 

k 2 - z 2 m% + m 2 _ 
Pi = ZiPij^ _ n . _ Pk + k±, (59) 

ft ~\~ i kj-(l-^) 2 m?-+m| _ 
* = + (1 - ~ Zl ) 2 PlJ p k Pk ~ k ^ (60) 

;i -ti -/^ 2 )( i -Vij,k) 



Ffe = ( 1 =71 Pfc , ( 61 ) 



with the spacelike transverse-momentum vector k± pointing in a direction perpendicular 
to both the emitter and the spectator momentum. 
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3.1.2 Massless case 

The case of a final-final splitting is considerably simpler in the massless limit, i.e. where all 
occurring partons can be treated as massless, — p 2 , — p 2 . = p\ = p 2 = 0. In this case, of 
course, the variables chosen to specify the splitting remain unchanged with respect to the fully 
massive case. However, neglecting masses the ordering parameter reduces to 

k l = Q 2 Vij,k Zi(l - Zi) = Ipijfik Vij,k Zi (1 - Zi) , (62) 

with the identification of Q 2 = 2pijpk this is identical with the transverse momentum defined 
in Eq. ([9]). The full phase space for the emission of an extra parton extends to Z{ G [0,1], 
Vij,k £ [0, 1], whereas again uniformly covers the interval [0, 2ir]. 

In the massless limit also the spin averaged splitting kernels (V^) simplify considerably, 
namely to 

(V^^, y ijtk )) = C F J t—^—z (1 + Si) \ , (63) 

OWfoywO) = 2cA + - 1 _- -2 + ^(1-^)1, (64) 

^ 1 6% \ ^iyij,k ' yij,k ^iyij,k ) 

CVW(^)) = ^{1-25,(1-*)}. (65) 
When combining the factorised form of the (m + l)-parton phase space, 

d$ m+1 = d$ OT £ £ T^l — d~Zi ^ (1 - Vij,k) e(Zi (1 - Si)) e(y ij)k (l - y ij>k )) , (66) 
with the corresponding expression for the [m + l)-parton matrix element, 

\M m+1 \ 2 = \M m \ 2 £ £ 7— T73^ (V ii)fc (Ii, , (67) 
ij k^ij -/'■I'j • v ij 

the fully factorised form of the (m + l)-parton differential cross section is recovered 

da m+1 = dcr m £ ^^dZi 7T 7T TTipec ^ (yij,k) (Vij,k{Zii 2/ij',fe)) ■ (68) 

y fc^- //'.;./■• 2;-i 2tt A/y 
However, in this case, the Jacobian J(yij ; k) simply is given by 

J (Vijje) = 1 - Vij,k ■ (69) 
With the transverse momentum defined according to Eq. ( 1621) again the identity 

dyij,k _ dk^ 



Dij,k 



(70) 



is found. Choosing k 2 ± as the evolution variable with its lower cut-off given by k^ and the 
upper limit by kj_ max the £j integration range reduces to 



1 



"^(kj^max) ^Lfl) — -I 1 =F a J- , 2 

\ \ l,max 
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k 2 

K _i_,o 



(71) 



Figure 2: Sketch of the splitting of a final- state parton accompanied by an initial-state spectator. 
The blob denotes the m-parton matrix element. The incoming and outgoing lines label the 
initial- and final-state partons, respectively. 



Given a valid set of k 2 and Si this can be solved for 

k | . . 

Vll ' k = Q*zi(l - z{) ' (72) 

completing the determination of the splitting variables. Making the necessary replacements 
when going from massive partons to massless the Sudakov form factor given in Eq. (|48|) yields 
the corresponding non-branching probability. The massless kinematics can be derived from 
Eqs. (l59l-(l6TT) by setting /i^- = /ij = \ij = 0, accordingly 

k 2 

Pi = ZiPij-\ , ;jv , — pk + k±, (73) 

Zi 2pijp k 

» = (1 - i ' ) ^ + (rfc"'" tl ' (74) 

p k = (l-y ij)k )p k . (75) 



3.2 Final-state emitter and initial-state spectator 

In this section, the case of a final-state emission with the spectator being an initial-state parton a 
is worked out. The splitting schematically reads {ij, a} — > {i, j, a}, for a pictorial representation 
of the configuration, cf. Fig. [2j This configuration emerges for the first time when considering 
deep- inelastic lepton scattering (DIS), where one incoming line carries colour charge, or in 
configurations like vector boson fusion, with no colour exchange between the two hadrons. 
However, besides the singularity related to a final-state splitting, there is also a singular region 
for the splitting of the initial-state QCD parton, which needs to be included in such processes. 
This situation will be investigated in detail in Sec. 13.31 



3.2.1 Massive case 

The initial line is always assumed to be massless, however, all final-state particles can be 
massive. Accordingly, 

~2 2~2 2n2222 

Pij = m tj p a =p a = Q p. = m 8 , p j = mj . (76) 
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To avoid on-shell decays being described incorrectly, again m 2 ,- < mf + m 2 should hold true. 

• Kinematics: 

Four-momentum conservation is incorporated through the condition 

Pij -p a =Pi+ Pj ~Pa = Q ■ (77) 

Defining the Lorentz-invariants 

PiPa + PjPa - PiPj + limtj - m f - m 2) 

X%j a . > I ' o ) 

PiPa + PjPa 

z- = 2^1 z- = ^ = l-z- (79) 

PiPa + PjPa PiPa + PjPa 

the relative transverse momentum of the new emerging final-state partons is given by 

ki = 2p a p ij 1 ~ Xij ' a ~z t (1 - ~z t ) - (1 - z t ) 2 m 2 - 5> 2 . (80) 

Xij,a 

• Phase space: 

The factorised form of the three-parton phase space reads [3T] 

i 

d$(pi, Pj; Q+Pa) = J dxd^ipif, Q + xp a ) [dpitpijipa, x)] 0(x + 



with the single-parton phase space factor 

[dpiip^pa, x)) = 2 f* jP " ^ d^ clary a 5(x - x ija ) , (82) 
and the integration boundaries 

x- = , x + = 1 + //J - (/ii + /ij) 2 , (83) 



1 - x + ii% + nl - \i) T - x + [i% - /if - p 2 ) 2 - 4/i 2 /x 2 
% " 2(1 -x + 4) • (M) 

Here, again rescaled parton masses have been introduced, 

y-n = i m ™ = ( n = h j, ij) • (85) 

Sj ^PijPa/ Xij ja 

Splitting kernels: 

The polarisation-averaged QCD dipole splitting kernels (V^(zi, Xy^)) read 

f 2 m 2 1 

<V S(gJ (^,)> = Cf + , (86) 



^ij,aj PiPj 

+ r -2 + Zi(l-Zi 



%i (1 Xij <a ) Zj -|- (1 -£ij',a) 

(87) 

(V^ iQ .(^)) = T R {l-2(z+ -!,)(«_ -I,)}. (88) 

The scalar product of the a priori unknown momenta pi and Pj in Eq. ( !86l) can again be 
expressed according to Eq. ( HOI ). 
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Matrix element: 

Combining the (m + l)-parton phase space with the factorised form of the matrix element, 
\M m+1 1 2 = \M m \ 2 £ £ 7— —4— -2 T7^ — - 8^s (V^.(^, xy, „)> , (89) 

one obtains the fully differential cross section for the emission of one additional parton in 
that configuration 

da m+1 = da rn J2J2 7T" Si T7^ : * (Vy&, ^>)) , (90) 

where the sum covers all the possible colour-connected emitter-spectator pairings. The 
Jacobian of the variable transformation in this case reads 

Jfaj,a) = ^rr — • (9i) 

Taking into account that the initial parton actually stems from a hadronic initial state, 
a corresponding parton distribution function (PDF) emerges. Absorbing it into the Ja- 
cobian yields 

A^,a, /4) = i ~ l^hl^M . (92) 



x. 



ij,a 



Here, r\ a is the momentum fraction of the spectator parton a and f a (Va, Hf) the corre- 
sponding hadronic PDF evaluated at some scale [xp. In Sec. 12.31 this scale has been set 
to hf = kj_. The parton distribution function f a (Va/ %ij,a, Hf) corresponds to the new 
incoming momentum and is also evaluated at scale fip. 



Sudakov form factor: 

Note that Eq. Mjh implies that 



dxij ;a /i N dk 



2 



Xij,a 



(1 x ij,a) , 2 ' (93) 



With taken as the evolution scale with an upper limit k^ max and the cut-off k^ the 
Zi integration boundaries therefore are given by 



^-(kL^kin) = Maxlill 



z+(K^Ko) = Mini 1(1 + 



k 2 

K _i_,o 



2 

_L,max 



2 I \ k 



k 2 



,2 

_,max 



z \ , (94) 
(95) 



with z± given in Eq. ( !83l) . Having determined k^ and Zi the variable Xij A is calculated 
through 



2\ 



_ kj + (1 - Zifrnj + z\m) - ~z % {\ - Zj){m\ - mf - m 3 , { , 

%tha ki + (1 - -ZiYrnl + 5 2 m 2 + Zi{l - ~ Zi ){Q 2 + 2m 2 + 2m 2 ) ' 1 j 
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and has to fulfil the condition 

X%j,a e [Va/Vmax, X+] (97) 

to yield a valid branching. Here, ?7 max corresponds to the maximal allowed Bj0rken-x 
for the PDF. Having at hand all ingredients, the Sudakov form factor associated to the 
splitting of a final-state parton with an initial-state spectator reads 

AFi(k^ max , k^ )0 ) 

\ \^ 1 f dk ± r a s (k ± ) f a (r] a /xij ta ,k ± ) , XTa( ~ ^ 

J 



exp 



A^ ec J kl J 1 2vr f a ( Va ,k] 



k 2 „ 2 - 



(98) 

Physical kinematics: 

The actual branching kinematics can be given in a Sudakov parametrisation. In the 
Breit-frame of the emitter and spectator the two final-state momenta can be written as 



k 1 + nij — ;// 



~"1 ml 



Pi = ZiPij^ 9 - 1 Pa + k±, (99) 

with the sp acelike- k± being perpendicular to both the emitter and the spectator mo- 
mentum. After the splitting the latter remains parallel to p a but is rescaled according 
to 

Pa = —Pa- (101) 

3.2.2 Massless case 

The modifications emerging in the massless limit are briefly discussed. The splitting variable 
Xij ta simplifies to 

PiPa+P^-piP^ (1Q2) 
PiPa + PjPa 



whereas the momentum fractions Z{ and Zj are still defined according to Eq. ( 1791) . The invariant 
spacelike transverse momentum is simplified and reads 

ki = 2p a p ij 1 - Xij ' a ~Zi{l-~ Zi ). (103) 

While the g — > gg splitting function remains the same, the mass dependent terms drop out 
in the q — > qg and g — > qq kernels, 

(Vl gj (z h x ijia )) = (J F | — ^— ^-^ -(i + ^j, (104) 

(v; qj (5,)) = r R {i-2^(i-^)}! (los) 
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Incorporating the factorisation of the (m + l)-parton matrix element and the corresponding 
phase space the fully differential (m + l)-parton cross section is still given by Eq. f[9"01 . with 
the appropriate Jacobian for hadronic initial states. In the massless limit the phase-space 
boundaries are no longer constrained through finite mass terms, and therefore extend to 



■Kij,ai Z{ G [0, 1] . 



Eq. ( EEDHD still implies that 



dXjo a 



(1 •Eij,a) 



dk 2 ± 



(106) 



(107) 



When evolving in from max and asking for a minimum separation k^ the allowed Zj 
range is reduced to 

Zi G 



\ 



1 



ill 



\ 



in the massless case. The expression of the Sudakov from factor, Eq. 
unaltered. 

The kinematics of the new final-state partons simplify to 

k 2 

Pi &i Pij ~r 



Pi 



[1 - Zi)pij + 
b 

spectator momentum is still given by 



Z% ^PijPa 

T^^P 



Pa + k A 

k 



(108) 
of course remains 

(109) 
(110) 



[I- 2pijp a 

with k± still being perpendicular to both the emitter and the spectator momentum. The new 



1 



Pa 



X 



-Pa 



111) 



ij,a 



With Xij ja taken from Eq. ffT02l . 



3.3 Initial-state emitter and final-state spectator 

The case of an initial-state parton branching (ai), accompanied by a final-state spectator (k) is 
sketched in Fig. [HJ This accounts for the situation where the emitter and the spectator parton 
studied in Sec. 13.21 exchange their roles. 



3.3.1 Massive case 

As stated above, treating initial-state particles as massless, final-state particles emitted from the 
initial state are assumed massless as well, the spectator mass, however, is arbitrary. Accordingly, 
the momenta involved in the splitting {ai, k} {a, i, k) have to fulfil the mass-shell relations 

P 2 ai=P 2 i =Pa = 0, p\=p\ = m\. (112) 
and the momentum conservation condition 

Pk-Pai=Pi+Pk-Pa = Q ■ (113) 
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Figure 3: Splitting of an initial-state parton accompanied by a final-state spectator. The blob 
denotes the m-parton matrix element. The incoming and outgoing lines label the initial- and 
final-state partons, respectively. 

• Kinematics: 

The splitting can be specified by the variables 

PiPa + PkPa ~ PiPk PiPa 

Xika = , U{ — . (114; 

PiPa + PkPa PiPa + PkPa 

The transverse momentum squared parametrising the singular region where the emitted 
parton i becomes collinear with the initial-state parton a then reads 

k]_ = 2p ai p k - — Ui(l - Ui) . (115) 

ik,a 



fi k = , mk : (116) 



To allow for a more compact notation, the rescaled spectator mass 

\j2PaiPkl 'Xik,a 

is introduced. 
Splitting kernels: 

The QCD splitting kernels, taking into account possible non-zero spectator masses, read 

(vri^,^.)) = c F \ 2 — (i +Xikta )\ (H7) 

^ A %ik,a ~T Ui J 

(VTW) = C F L k , a + 2^^-^-^-\ , (118) 

^ik,a ■Eik,a J- Ui J 

(Vf g ^, a , Ui )> = 2C A { 1 +j_ ^fea-i 

J- %ik,a i" %ik,a 

+X ik , a (l ~ ***,«) ~ — ZT^-) , (H9) 

(V^^m)) = T K {l-2x iKa {l-x iKa )} . (120) 
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Phase space: 

The three-parton phase space is again obtained by a convolution of a two-parton piece 
and a single-parton part, 

i 

d®(pi,p k ;Q +p a ) = J dxd<&(p k ;Q + xp a ) [dpi(p k ; p a , x)) , (121) 



where 

1 X! 2 



[dpi(p fc ;p a ,x)] = — — 8{p i )Q(x)Q{l- x)5{x -x ika ) - , (122) 

Z7T 1 — Ui 

or, more conveniently, 

[dpi(p k ;p a ,x)] = ^^^dx ifcia dMie(Mi(l-Mi))e(x(l- x))5{x - . (123) 
The upper limit for the Wj-integration contains a dependence on the spectator mass, 



u+ = 1 Xtk ' a - . (124) 

Matrix element: 

Using the factorisation property of the (m + l)-parton matrix element 

\M m+ l\ 2 = \M m f Y,T,7^J^E^—^s(yt^ik,a,U l )) (125) 
ai k ^VaVi • /V oi x ik,a 



in the soft and collinear limits and the relation 

the (m + l)-parton fully differential cross section reads 



(126) 



da m+1 = da m — dx ^,a ^ TP TTipec ~ ^ (Vfc*(Xifc,a, ■ ( 127 ) 

The integration range of the variables Ui and Xi k<a is [0, u+] and [0, 1], respectively, and 
[0, 27r] for 0. The Jacobian 

J{x lk>a ) = — (128) 

for the parton matrix element again is changed in hadronic interactions to include the 
effect of the PDFs, such that 

J(x lk , a , ft F ) = — , (129) 

%ik,a Jaiylaii ftp) 

where again, in the implementation here the choice for the factorisation scale is \xp = kj_, 
cf. Sec. 12.31 Note that the Jacobian takes into account not only a change in Bj0rken-a; 
but also a possible flavour change in the process' initial state. 
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Sudakov form factor: 

The integration over it, in Eq. (\V27\i can be replaced by an integration over according 



to 



Ui 



l — Ui dk 2 ± 
1 - 2 Ui 



(130) 



The arising Jacobian is combined with the function J(xik, a ', a40 to 

J{X' U ' jJ^ ) ^ ^ faijlai/ -^ik,aj ^f) 

1 2ltj 3^ifc,a fai{Vaii A*f) 

With > as the evolution variable and its cut-off being kj_ the a;^ phase-space 
boundaries are 



(131) 



•^ik.a ^ 



'at 



//max 4" 4k 2 L 



(132) 



with r/max the maximal allowed Bj0rken-x of the PDF. With k[ and Xik, a given, u% can be 
calculated and yields 



" = 2 " 



Ahj_Xik t a 

\" ~ g 2 (i - X ik , a ) 



(133) 



When Ui < u + an allowed branching is found. Thus the Sudakov form factor for having 
no emission from an initial-state parton accompanied by a final-state spectator between 
scales k 2 max and k^ can be written down, 



Aip(k 2 i „„.k 



,max' _L,0/ 



_L ,max 



exp 



EE 



Kfspec 
ai k yV tii 



dk 2 ^ 



dXjfc, a - J{%ik,a-i u i] kj_) (V fc (Xik,ai u i)) 



k 2 

K ±,0 



/ 

(134) 



• Physical kinematics: 

The new initial-state particle a remains parallel to the original initial-state parton, and 
is just rescaled by the splitting variable Xik, a such that 

Pa = -^—Pai- (135) 

ik,a 

The two final-state momenta are most conveniently evaluated in the rest-frame of Q + p a 
with p a pointing along the positive z-axis. The corresponding light-cone momenta read 

Q + Pa = {M,M,6) and p a = {2E a , 0, 0) . (136) 



26 



Note that the massless vector p a only has a light-cone +- component, given by twice the 
energy of the parton. For pi and pk the ansatz 

— * — * 

Pi = (l±e y ,l ± e~ y ,l ± ) , p k = (m k:± e z ,m k:± e~ z , -l±) , (137) 

is used, with m± being the transverse mass. Besides the energy- and momentum- 
conservation requirement the momenta are constrained by the splitting variables, 

PiPa l±e~ y , . 

Ui = ~, : = — n— ■ 138 

yielding 

P ± = (M 2 - m 2 k )u t - M 2 u\ , (139) 

for the transverse momentum squared. This equals the physical transverse momentum 
squared of parton i, k^. Employing the relations 

M 2 - ml . 1 f l± Mui \ 

coshy = -mr* sinhy = 2{mrr~L) ' (140) 

cosh, = fj m «, sinh^lf - M ^-^\, ( 14 1) 

2Mm k:± ' 2^M(1-^) m fc)± ^ ' 1 j 

the four-momenta of the final-state partons, in the frame specified above, read 

Pi = (l± cosh y, l± cos 0, l± sin 0, l± sinhy) , (142) 
p k = (rrik ; ± cosh z, —l± cos 0, — Zj_ sin 0, m^j. sinh 2) . (143) 

Again, has been uniformly distributed in the transverse plane. The kinematics is 
completed by rotating and boosting the momenta p a , pi and pk back in the laboratory 
frame. 



3.3.2 Massless case 

The massless limit of the scenario above, initial-state splittings accompanied by final-state 
spectators, {ai, k} — > {a, i, k}, corresponds to neglecting the spectator mass, p\ — p\ — 0. 
Apart from that, the splitting variables remain unchanged and the dependence on m&, of 
course, disappears in the corresponding phase space boundaries. 

Dropping the explicit mass terms present in (V^° qi (xjfc ia )) and (Vf agi (x;fc ia , ui)) given in 
Eqs. ( [TTBJ) and (1119p . respectively, the factorised form of the fully differential cross section can 
completely be taken over. 

Neglecting the finite spectator masses the splitting kinematics is significantly simplified. In 
the emitter-spectator Breit-frame 

Pa = —Pai, (144) 

%ik,a 

Pi = (l-Ui)- — ^^p ai + Uipk + k ± , (145) 

%ik,a 

Pk = Ui - — X ^ p ai + (l- Ui) p k - k± , (146) 

ik,a 

with k± perpendicular to both the emitter and the spectator. 
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b / / p h 



Figure 4: Schematical view of the splitting of an initial-state parton with an initial-state parton 
as spectator. The blob denotes the m-parton matrix element. Incoming and outgoing lines label 
the initial- and final-state partons, respectively. 

3.4 Initial-state emitter and initial-state spectator 

The last scenario to be studied is the splitting of an initial-state particle ai, with the spectator b 
being an initial-state parton as well, cf. Fig. HI This type of branching occurs when considering 
hadron-hadron collisions, where both the initial-state particles are colour charged and therefore 
can be colour connected. The simplest example for this configuration is the lowest order Drell- 
Yan process, where both the incoming quark and anti-quark can serve as emitter and spectator. 

In contrast to all other cases discussed before, it turns out to be convenient to preserve the 
spectator momentum p in this branching. Since also the emitter momentum remains parallel 

tO p a , 



the transverse momentum of the emitted parton, p^, has to be balanced by another final-state 
momenta kj. This does not only include the QCD partons, but all non-QCD particles, e.g. 
leptons, as well. 

• Kinematics: 

Defining the variable 



with Xi, ab 



PaPb ~ PiPa ~ PiPb 



(147) 




PaPb 



PiPa 




PaPb 




(149) 



The four-momenta of the m-parton ensemble fulfil 



in 



Pai + Pb-J2 k 3 =0 i 

correspondingly the full set of m + 1 particles has to satisfy 



(150) 



711 



Pa + Pb ~ k i ~ Pi = ' 



(151) 
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Splitting kernels: 

The polarisation-averaged splitting kernels (V aj ' b ) depend on x^ only and read 

2 'l + x iiab )\, (152) 

(153) 

l + x i>ab (l -x itab ) \ , (154) 





( 

= c F l 


(V q ° q " 6 (x^)> 


( 

= c F \ 


(V^ix^)) 


= 2C A 


(V^ b (x i;ab )) 


= T R { 



1 Xi,ab 

%i,ab ) 
I 1 ,1 ^i,ab 



1 ^i,ab ^i,ab 



(155) 



Phase space: 

The final-state phase space can be written as follows [31] 

i 

d®(p h k 1 ,...;p a + p b ) = J dxd§(k 1 ,...;xp a +p b )[dpi(p a ,p b ,x)}, (156) 

o 

with 

[dpi(p a ,p b , x)] = ^-dx i:ab dviQ(x(l - x)) Q(vi) ( 1 - r-^—) S(x - Xi :ab ) , (157) 

l07T z Z7T \ 1 — XJ 

where is the polar angle in the plane perpendicular to p a and p b . 
Matrix element: 

Combining this with the expression for the (m + l)-parton matrix element 

\M m+l \ 2 = \M m \ 2 £ 7^— -tfsp^ — 87ra s (V ai ' b (xi,ab)) (158) 

ai b^ai ^VaPi ai X ljab 

the differential cross section becomes 

da m+1 = da m dx ^ 7T 7T 77^ ~~ ~~ ( vahb ( x i,ah)) , (159) 

ai b+ai V i 271 271 N ai X i,ab 

where 1 — Xi >ab — V{ > has to hold. The Jacobian can be read off as 

J(Xi,ab) = — , (160) 
•Ei,ab 

or, including again the PDFs, 

J (**,«&; ^f) = 7 7- 2 x • (161) 

•^ijOfe J aiy'laii H'F) 
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(162) 



Sudakov form factor: 

Regarding the transverse momentum given by Eq. (j!49p the identity 

dvj _ 1 - x i<ab - Vj dkj_ 
Vi 1 - x i<ab - 2vi k]_ 

can be employed to replace the Vj integration with a k^-integral. The resulting Jacobian, 
combined with J(x ija b; fi 2 F ), amounts to 

1 — £i,ab V{ 1 fa{Vai/ Xi,ab} l^p) 



1 Xi.ab 2fj %i,ab fai[?]aii I^f) 



(163) 



When evolving in the dependence of the x^-integration boundaries on the cut-off 



kj_ q read 



%i,ab ^ 



ZPaPb 



^max ' tyaPb + 4k^_ „ 



Vi can be calculated from k^ and x^, 



1 - Xi 



i,ab 



Vi 



The Sudakov form factor then reads 

,2 1,2 



2kU 



i,ab 



\ PaPb{^ ~ Xi, 



ab) 



(164) 



(165) 



An(ki imax ,ki i0 ) 
( 

= exp 



EE 



v 



Kfspec 
ai b^ai J v ai 



dk^ 



d.r, Jlb Qs( ^ /4) J(X <)B6 , T5 i; ki) (V^Zi.ab)) 



2tt 



(166) 



Physical kinematics: 

The momenta of the (m+l)-parton ensemble, expressed through the emitter and spectator 
momentum and the momenta of all other final-state particles of the m-parton process, 
read 

1 



Pa 



Pi 



Pai 5 

Xi t ab 

1 Xi ab Vi „ _ 
Pai +ViP b + k j 



fcj = A(Pai + Pb, Pa+Pb- Pi) kj , 

with \f^_ uniformly distributed in the transverse plane and A(p ai + p b ,p a + Pb~P 
A(K, K) being a proper Lorentz transformation given by 



(167) 

(168) 
(169) 



^ _ 2{K + KY{K + K) U 2K»K U 

+ k 2 



A» V (K, K) = g\ - ! ^ + . (170) 

(K + K) 2 

Accordingly, the full set of final-state momenta compensates for the transverse momentum 
of Pi, although they do not participate in the splitting. 
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Figure 5: The SUSY QCD vertices corresponding to gluon emission off (anti-)s quarks and 
gluinos. 



In the minimal supersymmetric extension of the Standard Model the sector of strongly inter- 
acting particles is extended by the superpartners of the ordinary quark- and gluon-fields [61J. 
The new particles participating in the strong interaction are the scalar-quarks, called squarks 
and the gluino. While the former are colour-triplets the gluino is a Majorana fermion in the 
adjoint representation, a colour-octet. 

In order to be consistent with todays experimental (non-)observations the assumed SUSY 
particles have to be rather heavy. This renders the massless limit for these fields not applicable 
when describing their QCD interactions at the energies of the forthcoming colliders. Based 
on that argument it is beyond the present scope to describe possible branchings like g —>■ qq*, 
g — > gg in a quasi-collinear limit. Rather, they are appropriately described using exact matrix 
element methods, as discussed e.g. in [681(69] . 

Since the spin and the flavour of the spectator parton do not enter the splitting functions, 
the branchings of the Standard Model particles are not altered in supersymmetric extensions. 
The only SUSY QCD splittings that appear to be relevant in the context of a parton shower 
formulation are related to the emission of a gluon off a squark or anti-squark and off a gluino, 
cf. Fig. El Further, assuming that supersymmetric particles do not appear as partonic initial 
states those are solely final-state splittings. The associated spectator, however, can be either 
in the final state or in the initial state. 

Due to its fermionic nature the splitting functions involving gluinos are equal to the corre- 
sponding splittings of massive quarks, cf. Eq. (j57|) and Eq. flHBT) . only the colour factors have 
to be adopted from Cp to Ca- 

The kernel of the branching q — > qg with the spectator also in the final state reads 



where all the variables have been defined in Sec. 13.1.11 If the spectator is in the initial state 
this becomes 



3.5 SUSY QCD splitting functions 



1 Z% ~\~ ZiDij,k ^ij,k 




(171) 




(172) 



for the definitions of the variables 



used see Sec. 13.2.11 
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Apart from the splitting kernels all the results derived in the corresponding sections describ- 
ing the branchings of massive final-state partons with spectators in the final- or initial state 
can be taken over without any alteration. This includes the exact phase-space factorisation as 
well as the parton kinematics defined there. 



4 Comparing the hardest emission with matrix elements 

In the following, the predictions for the hardest (first) emission of the parton shower algorithm 
will be worked out for different processes and compared with corresponding exact tree-level 
matrix element calculations. The set of processes to be considered covers three-jet production 
in e + e~ collisions, cf. Sec. 14.11 the first order real correction process to DIS, cf. Sec. 14.21 and the 
production of a weak gauge boson accompanied by a light jet at hadron colliders, cf. Sec. 14.31 
These three examples constitute a full set of generic processes to reliably test the first emission 
of the proposed parton shower approach. 



4.1 Three-jet production at lepton-colliders 

In this example the production of three jets at a lepton-collider is investigated. Jet production 
proceeds via the s-channel exchange of a colour- singlet particle, namely a 7* or Z°-boson. The 
latter will be ignored in the discussion here. At first perturbative order in a s , two Feynman 
diagrams contribute to the matrix element 7* — > qqg, corresponding to the emission of a gluon 
from either the final-state quark or the anti-quark, cf. Fig. [HJ 



For convenience, the centre-of mass energy 

E c . m . = ^, (173) 

and the momentum fractions 

(174) 

are introduced. Neglecting the masses of the final-state particles the Lorentz-invariant Man- 
delstam variables for the 1 — ^ 3 process become 

s = ( Pl +p 3 ) 2 = 2p lP3 = Q 2 (l-x 2 ), (175) 



32 



t = (p 2 +p 3 ) 2 = 2p 2 p 3 = Q 2 (l-x 1 ), 
u = (pi+ p 2 ) 2 = 1p\P2 = Q 2 {1- £3) • 

Energy-momentum conservation implies that 

x\ + %2 + X3 = 2 and s + i + u = Q 2 . 



(176) 
(177) 



(178) 



The partonic differential cross section with respect to the quark and anti-quark momentum 
fractions xx,i reads 



da 



dxidx 2 



ME 



Z7T 



^1 - X X )(1 - X 2 ) 



where &o denotes the total cross section of the two-jet process 7* 



qq, 



(To = 2a qed e 2 E, 



qJ-^c.m. ; 



(179) 



(180) 



see for instance [70 J. 

In the parton shower approach, two contributions occur as well. They correspond to the 
timelike splitting of either the quark or the anti-quark, and the total result is just the incoherent 
sum of the two pieces. To work this out, consider the case of the quark splitting with the anti- 
quark being the spectator parton. Then, the shower variables are, cf. Sec. 13.1.21 



2/13,2 



P1P3 



P1P3 + P1P2 + P2P2, 
P1P2 u 



s + u + t 



s 

Q 2 



P1P2 + P3P2 u + t 
which, expressed in terms of the x,, translate into 

1 — X3 

2/13,2 = 1 - x 2 and zi -- 



1 — Xi 



x 2 x 2 

Accordingly, the cross section for the emission off the quark can be cast into the form 



da 



dxidx 2 



PSq 



Co — o F 

Z7T 



1 



1 — X2 \2 — x\ — X2 



1 — X\ 



X 2 



181) 
182) 

;i83) 

184) 



The result for the emission of a gluon off the anti-quark can be obtained from Eq. f)184p by 
1^2. Taken together, the total parton shower cross section yields 



da 



dx\dx2 



da 



PS 



dxidx 2 



+ 



da 



PSq 



dxidx 2 



Z7T 



/y>2 I ,-vt 2 



'1 - Xi)(l - x 2/ 



+ 



PSq 

1 — Xi 

X2 



1 - X 2 

X\ 



185) 



Obviously, the parton shower cross section reproduces the matrix element calculation in both 
the soft and the collinear limit. The only difference between the two results are two non-singular 
terms in the parton shower result that vanish as x\^ — > 1. 
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7*/Z (g) 



q(pi 





q(p) g(p2) q(p) g(p2) 

Figure 7: T7ie two leading order Feynman diagrams contributing to 7*q — > 



tVz°(?) 



q(pi) 



tVz°(?) 



qfe) g(p) q(pi' 

Figure 8: T/ie two possible Feynman diagrams for 7*g — > qq. 



4.2 Real corrections to leading order DIS 

The simplest physical process involving initial-state hadrons is deep-inelastic lepton-nucleon 



scattering (DIS), i.e. e^p 
namely e ± q — > e ± q and e ± q 



e + X. At leading order, two partonic processes contribute, 
♦ e q, both of which must be convoluted with the initial hadron's 
PDF to obtain the hadronic cross section. The interaction is mediated by virtual-photon and 
Z°-boson exchange. In the following, however, only the 7* channel is taken into account, for 
which the two partonic cross sections are equal. 

At next-to-leading order the quark can radiate a gluon before or after its interaction with 
the virtual photon, cf. Fig. [71 Beyond this, at NLO the incoming quark may originate from a 
gluon in the initial hadron that produces a quark-anti-quark pair which the 7* then couples to, 
cf. Fig. [HJ The real emission matrix elements can be expressed through the kinematic variables 



Q 1 



x 



2pq' 



PiP 
pq 



(186) 



where q denotes the four-momentum of the off-shell photon, p the incoming parton momentum 
and the pi label the momenta of the final-state partons. The Mandelstam variables for the 
2^2 processes 7*(<?)q(p) — > q(pi)g(P2) and 7*(?)g(p) — > q(pi)q(P2) are 



8 = (q + p) 2 
t = ( Pl -q) 2 



2pq - Q 2 
-2 Pl q - Q 2 



-Q 



X 



X 

2 X + Z\ 

X 



187) 
'188) 
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2 x + z 2 



X 



u = (p 2 - q) 2 = -2p 2 q -Q 2 = -Q 

Momentum conservation implies that q + p = p\ + p 2 and 

s + i + u + Q 2 = 0. 
In the following, the two real emission processes will be discussed in detail. 

4.2.1 The gluon emission process 

The matrix element of the gluon emission channel 7*(q , )q(p) — > q(pi)g(p2) reads [30|7U] 

x 2 + z\ 



(189) 



(190) 



\M 2 ^p Xl P2\p)\ 



2 
ME 



(l-a;)(l-zi) 



+ 2 l-3x^i 



\Mia(q + xp;xp)[ 



(191) 



with A4i tq (q + p;p) the matrix element of the lowest order process. 

In the parton shower approach two contributions to this final state emerge. First, the 
emission of the gluon from the initial-state quark with the final-state parton serving as spectator 
(IF) has to be considered. Second, the initial-state parton acts as the spectator of the final-state 
splitting q — >• qg (FI). 

• IF: 

The "parton shower" -matrix element of the initial-state splitting with final-state spectator 
is obtained from Eq. f)125p and is given by 



|- A ^2,q(Pl,P2;p)|p S if 
1 1 



2pp 2 X 2 l,p 



1 — £21,p + u 2 



(1 + £21,; 



\M hq (q + xp;xp)\ 2 , (192) 



where the appropriate splitting function, Eq. (11171) with ^1 = 0, has been inserted. 
Employing the identities 



£21, 



PlP + P2P - P2Pl U + t + § 



"2 



Pip + P2P 
p 2 p t 



1 



P2P + P1P u + t 

1 1 



u + t s + Q 2 
z 2 = 1 - Zx , 



2pp 2 x Q 2 pp2/pq Q 2 (l-z 1 ] 
the expression above becomes 



(193) 
(194) 
(195) 



|-M 2 ,q(Pl,P2;P)|p S if 



Q 2 (l- Zl ] 



2 — x — z\ 



[1 + x] 



IM^q + xp-xp)] 2 .(196) 
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FI: 

In full analogy the shower expression for the final-state emission process yields 



|A*2,q(pi,P2;p)|p Sfi 
1 1 



2ptp 2 x 12 , p 



1 - Zx + (1 - Xi2, p ) 



With 



X 12,l 



PlP + P2P ~ P\P2 



PlP + P2P 

this can be cast into the form 

,2 87ra s 



x and z\ 



PiP 



\Mi, q (q + xp; xp)\ 2 . 



PlP + P2P 



(197) 



\M 2 , q (puP2;p)\ 



PSfi 



Q 2 (l-x) 



C ¥ 



. 2 — x — z\ 



-a + x) 



\M liq (q + xp-xp)\ 2 , (198) 



where in addition 



2piP2 = Q' 



1 — x 



X 



(199) 



has been employed. 

Combining the two parton shower contributions yields the final result, namely 



|A^2,qbi,P2;p)|p s = \M 2 , q (pi,P2;p)\l Si{ + \M 2 ,q{pi,p 2 ;p)\ 



PSfi 



Q 2 



x 2 + z\ 



(l-aO(l-zi) 



• \M 1<q (q + xp; xp)\' 



(200) 



When comparing this with the exact perturbative result of Eq. (I19ip . it can be inferred that 
the parton shower exactly reproduces the soft and collinear singular structure of the matrix 
element as z\ — > 1 or x — > 1. The only difference is an additional finite non-singular term 
present in the exact result. 

4.2.2 The initial-state gluon channel 

Expressed in terms of the leading order matrix element the exact real emission next-to-leading 
order result for the process 7*(<?)g(p) — > q(pi)q(p2) reads [3111170] 



\M 2ig (p!,p 2 ;p)\ 



ME 



87ra s 



To 



(^ + (l-^) 2 )(x 2 + (l-a:) 2 ) 
zAl - z x ) 



+ 8x(l 



X 



\M^ q {q + xp;xp)\ 2 . (201) 



Starting from the leading order matrix element 7*(q , )q(p) — > q(pi) there is only one possibil- 
ity in the parton shower to reach the 2 — > 2 process, the splitting of an initial-state gluon 
into qq and the q interacting with the off-shell photon. The second matrix element diagram, 
corresponding to the interaction of the anti-quark with the 7*, here has no parton shower coun- 
terpart. However, when starting the shower from the charge conjugated leading order process, 
namely 7*(q , )q(p) — > q(pi), this contribution will occur while the 7*q interaction will be missing 
instead. The two terms are evaluated separately and then added incoherently. 
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Emission off the quark: 

The case of an internal quark propagator is discussed first. According to Eqs. (j!25p and 
( 11201) the parton shower approximation to the matrix element reads 
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Emission off the anti-quark: 

Starting instead the shower from the q initiated process, and emitting the quark into the 
final state yields, correspondingly, 
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Due to the charge conjugation invariance of the leading order matrix element, 

\M m {q + xp;xp)\ 2 = \Mi, q (q + xp; xp)\ 2 , (204) 
the two parton shower contributions can directly be combined and yield 

\M 2 , g {pi,p 2 ;p)%s = \M 2 , g (pi,p 2 ;p)\l S(l +\M 2 , g (pi,p 2 ;p)\l Sn 
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Again the parton shower matches the soft and collinear behaviour of the matrix element given 
in Eq. (12011) and reproduces the exact result up to non-singular terms. 

4.3 Associated production of a weak gauge boson and a light parton 

The lowest order production process of weak gauge bosons {W^, Z°, 7*) at a hadron collider 
proceeds via the s-channel fusion of two initial-state quarks. Without loosing generality W 
boson production will be investigated in the following. The leading order process then simply 
reads qq' — ► W^. At order a s there are three processes emerging: qq' — > W ± g, gq' — > and 
qg — > W ± q / . Considering on-shell W 1 * 1 bosons for simplicity^, only 2 — > 2 processes have to be 
discussed, which can be described using the Mandelstam variables 

3 = { Pl +p 2 ) 2 = 2p lP2 , (206) 
t = ( Pl -p 3 ) 2 = -2p lP3 , (207) 
u = (pa - p 3 ) 2 = -2p 2 p 3 . (208) 



5 This corresponds to neglecting the off-shell gauge boson decays which, however, do not affect the QCD 
dynamics of the processes under consideration. The decay products of the gauge boson can be introduced into 
the process using the narrow-width-approximation, or by incorporating the full off-shell propagator. 
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Figure 9: T/ie leading order Feynman diagrams contributing to the process qq' — > W ± g. 



Momentum conservation then implies that 

s + t + u 

where mw denotes the W ± -boson mass. 
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4.3.1 The gluon emission channel 

The first channel to be discussed is the gluon emission process qq' — > W ± g. At tree-level, 
there are two Feynman diagrams contributing to the matrix element, cf. Fig. [HI The partonic 
differential cross section can be written as [TDj 
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with &o the cross section of the leading order process qq' — > W = 

3 S 27T 



(211) 



In the parton shower approach there are two ways to produce the final-state gluon, which have 
to be added incoherently: either the gluon can be emitted from the initial-state quark or from 
the anti-quark. 

• Emission off the quark: 

As a first step, the kinematical variables used in the parton shower approximation should 
be related to the Mandelstam variables. In the first case, the initial quark as emitter and 
the initial anti-quark as spectator, the parton shower variables become, cf. Sec. E2 
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Using the appropriate splitting function of Eq. fll52p . the parton shower differential cross 
section 
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can be cast into 
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and multiplying with s/s yields 
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Emission off the anti-quark: 

Swapping the role of the emitter and the spectator parton amounts to only interchanging 
t and u in the results above. Accordingly, the differential cross section in this case is given 
by 
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The full parton shower result is the sum of the two contributions and reads 
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Again, the parton shower approach provides the correct description for soft and collinear phase- 
space configurations but misses non-singular terms. The difference of the parton shower and 
the exact result can be quantified by the ratio 
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which can take values between 0.5 and 1 in full agreement with the result of the parton shower 
algorithm implemented in Pythia [7TJ. This indicates that the parton shower approximation 
tends to overestimate the matrix element - a feature already present, e.g. , in e + e~ — ► qqg. 



4.3.2 The initial-state gluon case 

There are two Feynman diagrams, cf. Fig. [TO, contributing to the channel with an initial-state 
gluon, i.e. to the process gq' — > W^q. The result of the full matrix element calculation reads 
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Figure 10: TTie leading order Feynman diagrams contributing to the process gq' — > W^q. 



In the parton shower approach only one emission process contributes to this channel, corre- 
sponding to the t-channel diagram. The s-channel contribution is not realised in the shower 
ansatz. Using the definitions of the shower variables as given above and the corresponding 
splitting function, cf. Eq. (I155p . the parton shower cross section reads 



da 

di 



PS 



s 2vr R -t 1 



where 



1 



— 7T~ 1 n 

S 2,71 



t + u 



22:3,12(1 - 3:3,12)] 

s 2 + 2m^(i+u) 
-si 



(221) 



and x 3i i2 = rn^/s 
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has been used. The ratio of the parton shower and the matrix element result is 
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varying between 1 and 3 [7T]. Accordingly, the parton shower ansatz tends to undershoot the 
exact matrix element. However, the shower is constructed to give the correct answer in the 
logarithmically enhanced phase-space regions and thus has the correct limiting behaviour in 
the soft and collinear limits. The differences identified here are a result of differences in the 
non-singular terms, contributing only in hard regions of phase space. The process qg — > W ± q / 
closely follows the above example solely t and u have to be exchanged. This leads to the same 
qualitative results and the same conclusions. 



5 Applications 

In this section, the abilities of the newly developed parton shower formulation in describing 
QCD dynamics will be highlighted by comparing its results for various physics processes with 
experimental data and other calculations: In Sec. 15.11 the predictions for hadron production in 
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e + e~ collisions as measured at LEP will be studied and some results related to a future machine 
operated at y/s = 500 GeV will be discussed. In Sec. 15.21 emphasis is put on the capabilities 
of the shower to describe particle production at hadron colliders such as the Tevatron or the 
upcoming LHC. 

5.1 Jet production at e + e _ colliders 

Measurements of hadronic final states produced in e + e~ collisions provide a very precise probe of 
QCD dynamics in the final state and an excellent means to deduce its fundamental parameters 
such as the value of as(m^) , see for instance [72] , and the colour charges Cp and Ca in three- and 
four-jet events as discussed e.g. in [72]-[75]. Therefore it is not surprising that in the past years 
calculations for relevant three-jet observables, such as thrust, have become available at NNLO 
[76] and that full parton-level Monte Carlo codes for four-jet final states at NLO have been 
constructed [77778] . Obviously such observables also provide a critical test of the corresponding 
final-state radiation piece of a parton shower model. However, due to the fragmentation of 
partons into hadrons, which at the moment can be simulated with phenomenological models 
only, the parton shower predictions can not directly be compared with experimental data but 
rather have to be supplemented with a hadronisation model. The new parton shower presented 
here therefore has been interfaced to the Lund string fragmentation routines of PYTHIA version 

6.2 [79J in the framework of the SHERPA event generator. The large number of very precisely 
measured observables at various energies allows tuning the intrinsic parameters of the parton 
shower in conjunction with the phenomenological parameters of the fragmentation model. Such 
a procedure has been performed, for instance, for the new parton shower and fragmentation 
code in HERWIG++ [EH]- In principle, such a tuning is a very time-consuming and delicate 
procedure, see for instance [ST], deserving a publication in its own right. Recent developments 
to automatise the task of generator tuning and validation to a large extend are reported in 
[82J. Here, only a very limited tuning based on few parameters and observables only has been 
performed. The results of this tuning are presented in Sec. 15.1.11 In Sec. 15.1.21 the focus is on 
heavy-quark production at LEP1 and ILC energies to validate the treatment of finite parton 
masses in the shower model. 

5.1.1 Comparison with LEP1 data 

The most extensive data set available to validate QCD Monte Carlo predictions are LEP mea- 
surements at the Z Q pole. A selection of event shape variables, multiplicity distributions, 
differential jet rates, four-jet angle measurements and various particle momentum distributions 
have been used to select values for the unconstrained, phenomenological parameters of the sim- 
ulation, namely the value of the strong coupling constant at mz, the infrared shower cut-off 
kj^o and the three Lund string hadronisation parameters a (PARP(41)), b (PARP(42)) and a q 
(PARP(21)). For the results presented in the following, they have been fixed to a s (mz) = 0.125, 
kj_,o = 0.63 GeV, a = 0.33, b = 0.75 GeV -2 , and a q = 0.358 GeV, respectively. This yields a 
mean charged multiplicity per event of (N c h) = 20.87 at y/s = mz, in good agreement with the 
experimentally found value of (N ch ) = 20.92 ± 0.24 [83]. 

Figures fTTl to [141 show some exemplary results obtained with the new shower implementation 
compared to DELPHI LEP1 data at \/s — 91.2 GeV [H3]. 



41 



1 -Thrust @ LEP1 



Aplanarity @ LEP1 




Major @ LEP1 

1 F 1 — 1 — 1 — H — 1 — 1 — r 



Minor @ LEP1 




E 



10 




10" 1 



S 0.2 = 
■§ 0.1 = 

^ 
15-0.1 = 

7-0.2 = 
o 



— DELPHI 96 

— CS show. + Py 6.2 had. 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

m 



Figure 11: The event-shape variables 1— Thrust (1 — T), Aplanarity (A), Major (M) and Minor 
(m) in comparison with DELPHI data \8c 



In Fig. [TT] the new algorithm, denoted as "CS shower" in the following, is compared with 
some event-shape measurements by DELPHI |83j. The distributions of thrust, thrust-major, 
thrust-minor and aplanarity are displayed. The lower panel of each plot contains the bin-wise 
ratio (MC-data)/data, and the yellow bands show the statistical plus systematic error of the 
respective measurements. All the observables are sensitive to the pattern of QCD radiation 
probing both soft and hard emissions off the shower initiating qq pair. The Monte Carlo 
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Figure 12: The and p^ ut observables measured with respect to the thrust axis compared to a 
Delphi measurement J8$j . 



predictions agree very well with the event-shape data. There is some slight excess at very low 
1 — T corresponding to two-jet like events. This region of phase space however is very sensitive 
to hadronisation corrections and therefore dominated by non-perturbative physics. The same 
reasoning holds for the major and minor distributions at low M or m. 

The transverse-momentum distribution within and out of the event plane defined by the 
thrust and thrust-major axes, (p™) and (p^ ut ), respectively, are presented in Fig. [121 While 
is quite well modeled by the Catani-Seymour shower, p^ ut is significantly underestimated for 
values above 1 GeV. This tendency, however, is observed in other QCD Monte Carlo simulations 
as well [83J. 

In Fig.[13]the predictions for the exclusive two-, three-, four- and five-jet rates in the Durham 
algorithm [M] as a function of the jet resolution y^ ham are compared with data taken by the 
DELPHI experiment [HS]- They all exhibit a sufficient agreement with data within the experi- 
mental uncertainty bands. For the four- and five-jet rate the shower seems to underestimate 
the region of y^ ha,m ~ 0.001, however, this region is also affected by hadronisation effects and 
a more sophisticated tuning may provide an even better agreement with data here. The depen- 
dence on the choice of hadronisation parameters is even more pronounced for jet resolutions 
smaller than 0.001 where the results for the new shower preferably lie on the upper side of the 
experimental uncertainty band. 

The last observables to be considered are jet angular distributions in events with four jets. 
These observables can not be expected to be too well described by a pure parton shower as 
they should probe spin correlations of the produced partons. Such correlations, however, are 
not taken into account in conventional showers but require full matrix element calculations 
(eventually combined with a parton shower) to be completely taken into account [iTfSo] . In 
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Figure 13: The n-jet rates R n for the Durham jet algorithm as a function of the jet-resolution 
parameter y^ Tham . Data taken from J8h^ . 



Fig. [T31 the predictions for the Bengtsson-Zerwas [87] and the Nachtmann-Reiter [88] angle are 
compared with DELPHI data |85] for events with four jets at a jet resolution y^ rhs,m = 0.008. 
Both results agree surprisingly well with data. A similar level of agreement is observed for the 
other two prominent four-jet angles, 0:34 and the Korner-Schierhofz-Willrodt angle, that are 
not shown here. 



44 




Figure 14: The Bengtsson-Zerwas and Nachtmann-Reiter four-jet angles compared with DELPHI 
data J83j . 

5.1.2 Jet rates in heavy-quark production 

The leading order of heavy-quark production at lepton colliders also proceeds through an in- 
termediate 7* or Z° in the s-channel. Since pair production of top-quarks was outside the 
kinematical reach of LEP, only the production of bottom-quarks is available at these energies 
to discuss the treatment of heavy quarks in the new parton shower algorithm. At a future 
international linear collider (ILC), operating at or around ^/s = 500 GeV, pair production of 
top-quarks will play a key-role in the physics programme. This is also true for the LHC where 
top-quarks will copiously be produced and constitute a major background in nearly all searches 
for new physics. Therefore, a correct description of the radiation pattern of heavy quarks will 
be of enormous importance. As already hinted at in Sec. 13.1.11 radiation off massive quarks is 
suppressed with respect to the case mq = 0, also known as "dead-cone" -effect [IB]. The impact 
is however rather small when considering 6-quark masses of 4.8 GeV at collider energies that 
are much larger. To illustrate the impact of the finite 6-quark mass in the shower approach the 
Durham two- and three-jet rates for 66-production at LEP1 are presented in the left panel of 
Fig. [151 There, results are shown for the fully massive case (i.e. the mass has fully been taken 
into account in the splitting kernels, the phase-space boundaries and the splitting kinematics) 
and for the massless case are depicted. As expected, in the massive case both R 2 and R 3 
are slightly enhanced at low values of y^ ham , corresponding to the suppressions of additional 
radiation that turns a two-jet event into three-jet and a three-jet into a four-jet event at the 
scale of the emission. 

In the right panel of Fig. [TS] the same observables are presented, but this time for the 
pair-production of 175 GeV top-quarks at a 500 GeV ILC. Obviously, the finite mass has to 
be taken into account in the description of QCD radiation off top-quarks, since the differences 
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Figure 15: The exclusive Durham two- and three-jet rates in inclusive bb production at LEP1 
(left panel) and ti production at a 500 GeV linear collider (right panel). The solid curves 
correspond to fully taking into account the quark masses in the parton shower simulation while 
for the dashed predictions the finite masses have been neglected. 



with respect to the massless case can exceed an order of magnitude for the two-jet rate. 



5.2 Particle production in hadron collisions 

With the advent of the LHC era, the description and simulation of particle production processes 
at hadron colliders gained even more relevance. Due to the colour-charged partonic initial states, 
every hard process at hadron colliders is accompanied by initial- and subsequent final-state 
radiation. In the following, only two examples shall be considered to highlight the performance 
of the new parton shower model in such situations. First, the inclusive production of Drell-Yan 
lepton pairs, the simple most process that features initial-state emitter - initial-state spectator 
dipoles, and, second, QCD jet production are discussed. For the latter, besides looking at some 
inclusive two-jet distributions, three-jet observables sensitive to the inclusion of QCD colour 
coherence are considered and qualitatively compared with data. 

For all the predictions presented below, the CTEQ6L set of PDFs |89j has been used, the 
strong coupling constant has been fixed to a s (mz) = 0.118 with its running taken at two- loop 
level, in accordance with the choice in the PDF, and the infrared cut-off of the shower is chosen 
to be kj^o = 2 GeV. Hadronisation of the partonic shower final states is again accomplished by 
an interface to the Lund string routines of PYTHIA 6.2 |79j . 
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lepton-pair rapidity in eV+X @ Tevatron Run I 



lepton-pair pseudo-rapidity in e*e +X @ Tevatron Run I 




Figure 16: The rapidity- (left panel) and pseudo-rapidity (right panel) distribution of e + e 
Drell-Yan pairs produced in pp collisions at yfs = 1800 GeV. 



5.2.1 Inclusive gauge boson production 

The production of electroweak gauge bosons, e.g. W ± and Z° bosons, and their subsequent 
decay into leptonic final states, is one of the most prominent processes at hadron colliders due 
to their clean signature. Although very interesting in their own right, their inclusive produc- 
tion, i.e. their production together with additional QCD jets, represents a serious background 
to many other interesting processes, like, e.g. the production and decay of top-quarks or SUSY 
particles. Therefore, many theoretical efforts have been undertaken to predict gauge boson 
production as precisely as possible, both at fixed order in the strong coupling, see for in- 
stance [90]-[94j, or focusing on the analytical resummation of large logarithms from soft gluon 
emissions, see for example |95j-[99j. An important ingredient in all cases have been parametri- 
sations of the PDFs and a good perturbative control over their scaling behaviour, which by now 
is known at the three-loop level [lOOj . In addition, in the past few years, Drell-Yan production 
formed the testbed for approaches aiming at the combination of tree-level matrix elements with 
parton shower Monte Carlos [17f2QII27lll01lll02] . Parton shower Monte Carlos thereby have to 
deliver the correct description for the bulk of the events where the bosons are accompanied by 
rather soft emissions only. 

In the following, Drell-Yan production of 7*/Z° at Tevatron Run I energies is considered 
with the bosons decaying into e + e~-pairs. They are constrained to fall into a mass-window 
of 66 GeV < M e + e - < 116 GeV. The predictions of the new shower algorithm will directly 
be compared to results obtained with the matrix element-parton shower merging approach as 
implemented in SHERPA. To this end, an inclusive sample combining matrix elements for no 
extra emission and one extra final-state QCD parton has been generated with SHERPA version 
1.0.10. In the figures this sample will be denoted by "SHERPA 1.0.10 CKKW (0+1 jet ME)". 
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Figure 17: The k± differential jet rates d\ and d 2 in e + e~ + X production at Tevatron Run I. 

The discussion of the results starts with the rapidity- and pseudo-rapidity distributions of 
the produced lepton-pair, see Fig. [THJ As the shape of the former is already described well 
at the leading order, i.e. without any radiation, there is hardly any difference visible for the 
two results. The gauge boson pseudo-rapidity distribution however, only emerges when there 
is some additional QCD radiation. The radiation pattern, and especially the hardest emission, 
determines this leptonic observable. The pure shower result is in excellent agreement with the 
merged result, which contains the exact tree- level matrix element for the first hard emission. 
However, the shower distribution is somewhat lower at central pseudo-rapidity and slightly 
exceeds the merged SHERPA result for the two maxima around r] e + e - « ±4. These differences 
can be traced back to the lack of sufficiently hard radiation in the shower, which is constrained 
from above though the default shower start scale for this process, namely the invariant mass 
squared of the initial dipoles, M e 2 +e _. Below that scale, however, the parton shower can be 
expected to deliver reliable results, and in order to fill the phase space above that scale, matrix 
element-parton shower merging techniques should be added. 

The smaller amount of hard radiation can be further quantified by looking at the differential 
jet rates d 1 and d 2 for the fcj_-jet algorithm |103| . displayed in Fig. [171 These observables deter- 
mine the scales where the first {d\) and second (d 2 ) additional parton gets resolved as a jet from 
the core process. The results for the Catani-Seymour based shower and the merged SHERPA 
sample agree well for small cluster scales but, as can be expected, the shower is significantly 
lower for values of > mz- 

The last observable to be considered is the transverse momentum distribution of the lepton- 
pair. This distribution has been measured with high precision by the Tevatron experiments. 
Like the Drell-Yan pseudo-rapidity it is very sensitive to both soft and hard radiation accom- 
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Figure 18: The p± distribution of e + e Drell-Yan pairs in comparison with data from CDF at 
the Tevatron, Run I (TU^. 

panying the produced boson. Fig. [18] contains a comparison of the prediction from the new 
shower model with a CDF measurement [104] [§. The agreement between data and simulation 
is quite good up to py's of approximately 80 GeV. The upper-right part of Fig. [TBI contains 
a blow-up of the low transverse-momentum region of pt < 20 GeV, this time, however, on a 
linear scale. There, the parton shower describes the turn-on of the distribution quite nice, the 
actual peak, however, is slightly higher and a bit broader than seen in data. To describe the 
very low transverse-momentum region a Gaussian-smeared intrinsic k± was introduced, with a 
mean of 0.52 GeV and a width of 0.8 GeV. A more detailed tuning of these values combined 
with the shower cut-off k_|_ may yield an even better description of the distribution's peak. 
Above 80 GeV the parton shower dies off very rapidly due to its phase space being constrained 
by the choice of the starting scale, k^ max = M e 2 +e _. For illustrative purposes a prediction has 
been added where the start scale has been enhanced to 4Mg +e _. While the results at low pr 
do not change significantly, the distribution continues in the tail, thereby following the experi- 
mental data. But, of course, with this choice of parton shower starting scale, there is a similar 
drop-off of the distribution at scales of around 4Mg +e _. However, since there is no guarantee 
that the parton shower kernels do perform well enough at large scales, i.e. outside the soft- 
and collinear phase-space regions, it seems to be overly optimistic to stretch its predictions to 
such high scales. Instead, the parton shower description should consistently be improved by 
incorporating exact higher-order corrections. 

5.2.2 Inclusive jet production 

The most obvious QCD production process to look for at hadron colliders is inclusive jet pro- 
duction. However, from a theoretical point of view this is quite a complicated process. Besides 
tree-level calculations for up-to six final-state jets, so far, there merely exist full next-to-leading 
order results up-to three-jet production [3S], |105| - [TU8"] . Despite of strong efforts, culminating 

6 A comparison of the merged SHERPA prediction with this data has already been presented in [17] . 
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Dijet invariant mass @ Tevatron Run I 




Figure 19: Dijet mass Mdijet measured by D0 at Tevatron Run I f!15\/ . 



in evaluating the complete set of necessary matrix elements [109|- [TT3] and in developing meth- 
ods to isolate the infrared divergences in the real correction part |114] a full NNLO calculation 
for inclusive jet production has not been finished yet. Also, from the point of view of the parton 
shower presented here, jet production at hadron colliders is rather involved. This is because 
the 2^2 hard process will contain all possible colour connections between initial-state and 
final-state partons. Hence, QCD jet production constitutes a severe test of the entire shower 
algorithm. The input parameters for the simulations have been chosen as specified above. The 
starting scale of the shower, however, is related to the transverse momentum of the 2 — > 2 core 
process' outgoing partons, namely k5_ max = p±j. 

The first thing to be looked at is a very inclusive quantity, the dijet invariant mass. This 
has been measured by D0 during Run I [115] . The jets considered there have been recon- 
structed using a jet-cone algorithm with a cone opening angle of R = 0.7 in the r) — <f> space 
and with jet transverse energies above 30 GeV. Dijet candidates have then been subjected to 
the requirement that both jets satisfy |t7j| < 1.0. Fig. [T9l exhibits the resulting dijet-mass dis- 
tribution starting at M di j et > 200 GeV. It is a very steeply falling spectrum spanning six orders 
of magnitude in the mass range under consideration. To compare with data the result of the 
(leading order) simulation has been normalised to the cross section observed in experiment. In 
fact, the prediction of the proposed shower algorithm then is in very good agreement with the 
data and almost everywhere exactly hits the weighted bin centers. 

Another interesting observable when studying dijet events is the azimuthal angle between 
the two highest-pr jets. If there is no additional QCD radiation the two jets have equal 
transverse momenta and they are oriented back-to-back. Thus, in this case, their azimuthal 
separation A^dijot = \<Pi — 02 1 equals tt. In the presence of merely soft radiation the azimuthal 
angles remain strongly correlated, the strength of the decorrelation rises with the presence of 
additional hard radiation. Therefore, the dijet decorrelation provides a testbed for soft- and 
hard QCD emissions without the necessity to reconstruct further jets. Fig. contains the 
results of a recent D0 measurement for cone jets found for R = 0.7 [116] . The data fall into 
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Figure 20: Azimuthal decorrelation in dijet events measured by D0 at Tevatron Run II 11161/ . 

different ranges of the leading-jet transverse momentum and are then multiplied with different 
constant prefactors in order to display them in one plot. In all cases, the second-leading jet 
was required to have a transverse momentum pt > 40 GeV and both jets are constrained to 
the central-rapidity region, \yj\ < 0.5. The data are overlayed with the respective predictions 
of the Catani-Seymour dipole shower approach. The simulation agrees very well with the data 
over the whole interval of A0dijet spanned by the experimental measurements. This is a very 
satisfying result as it proves that the proposed shower formulation not only correctly accounts 
for phase space regions related to soft and collinear radiation but also yields qualitatively and 
quantitatively correct estimates for rather hard emissions as well. Furthermore, since this ob- 
servable is quite sensitive to model-intrinsic scale choices such as the shower start scale and 
scales entering the running coupling constant and parton density functions, this agreement 
proves that the defaults have been chosen correctly. 

The last item to be discussed are observables in QCD jet production at hadron colliders 
that are known to be sensitive to the correct treatment of QCD soft colour coherence in the 
parton shower simulation. Colour- coherence effects have been widely studied for e + e~ collisions, 
for an early review see e.g. |117j . They manifest themselves in the fact that soft emissions 
are forbidden outside a certain angular cone around the emitting particle's direction, known 
as angular ordering [4^150] . To account for this in shower Monte Carlos the phase space for 
allowed emissions has to be properly constrained. Within the HERWIG Monte Carlo for instance 
this is realised by evolving the shower in terms of cone-opening angles. While the situation 
for pure final-state showers is quite clear, in hadronic collisions the situation is slightly more 
complicated due to the presence of more colour flows, among them those that connect initial- 
and final-state partons. As colour-coherence here already influences the first emission from the 
initial- and final-state partons QCD three-jet events are the best place to look for the pattern 
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of these phenomena at hadron colliders. 

In one of the pioneering studies |118] three-jet events that feature a hard leading jet and a 
rather soft third jet have been considered. Observables potentially sensitive to colour coherence 
are spatial correlations between the third jet and the leading ones. In |118j such discriminating 
variables have been found and by comparison with Monte Carlo simulations evidence for the 
observation of colour coherence in hadron collisions has been provided. This ultimately has 
led to a refinement of the PYTHIA shower algorithm in order to appropriately model colour 
coherence in the spirit of [119j . In the CDF study |118] jets have been defined through a cone 
algorithm with a cone radius of R = 0.7 and the following event selection criteria have been 
applied: 

• For the two leading jets the pseudo-rapidity is constrained to \r]i\ < 0.7 and \r]2\ < 0.7; 

• they have to be back-to-back within 20 degrees in the transverse plane, corresponding to 
|0i — 02 1 > 2.79 radian; 

• and the transverse energy of the leading jet, E T1 , has to exceed 110 GeV, the third jet is 
required to have E T3 > 10 GeV. 

• Only for the study of the a variable defined below the additional cut 1.1 < Ai? 2 3 < tt, 



A number of observables has been considered, the two most convenient and discriminating ones 
have been the pseudo-rapidity distribution of the third jet, and the polar angle in the space 
parametrised by A(p = 03 — 02 and AH = sign(^ 2 )('73 ~ namely a = arctan(Aif/| A0|) 0. 
It should be stressed that the published results, used for the comparison, are not corrected for 
detector effects, such as finite resolution and uninstrumented regions, and therefore can only 
qualitatively be compared with theoretical calculations. In contrast, the results of the Monte 
Carlo simulations exhibited in [118] have passed the full chain of the CDF detector simulation. 
In Fig. [2D the measurements are compared with simulated events at the hadron level. In the 
left panel the 773 distribution is shown and the right panel contains the comparison of the a 
distribution. Both predictions agree well with the data. The r] 3 distribution tends to be broader 
in models that take into account colour-coherence effects and only then theoretical calculations 
show the significant dip around r/3 « seen in data E. The a variable is also very sensitive to 
the inclusion of colour coherence. It decreases from a = —n/2 to a = but then the slope 
changes and the distribution rises as a —>■ tt/2. This trend is clearly seen for the simulation 
with the new shower algorithm. Models not taking into account coherence fail to describe the 
distribution's raise towards a — > tc/2 and have a clear excess of events at small \a\. Concerning 
the interpretation of these results the missing detector smearing for the shower simulation has 
to be kept in mind. However, in Ref. |118j estimates for the size of the detector effects are given, 

7 A further observable considered in the CDF study is the spatial separation of the second- and third jet in 
the r) — <f> space, Ai?23- This observable, however, seems to be less discriminatory between theoretical models. 
In addition, and more importantly, detector effects seem to have a larger impact on its discriminating power. 
Therefore it is not taken into account here. 

8 The HERWIG Monte Carlo, incorporating colour coherence through explicit angular ordering, describes the 
data very well. Switching on the approximate version of angular ordering in PYTHIA, realised by a veto on rising 
opening angles during shower evolution, significantly improves PYTHIA's agreement with data. 



where AR 
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Figure 21: The pseudo-rapidity distribution of the third- hardest jet (left panel) and the distri- 
bution of the angle a (right panel) in inclusive QCD three-jet production in comparison with 
CDF data taking during Tevatron Run I. Experimental errors are statistical only. Histograms 
are normalised to one. 



showing that the impact of the finite detector resolution is much smaller than the size of the 
physical effects. The generic features of the two observables presented here are not dependent 
on detector effects, and they are well described by the new shower formulation. 

The conclusion of this is that the proposed parton shower algorithm with its notion of 
emitter-spectator dipoles associated with the color flow of the event and using transverse mo- 
menta as evolution variable accounts for soft colour coherence and yields a very satisfying 
description, both on the qualitative and the quantitative level. It can be anticipated that such 
non-trivial quantum phenomena are of large importance at the LHC, since the phase space for 
jet production is much larger and hard jets are produced copiously. For a solid description of 
QCD therefore the systematic and correct inclusion of these effects is paramount. 



6 Conclusions and outlook 

In this publication a parton shower model based on Catani-Seymour dipole subtraction kernels 
has been presented, which was proposed for the first time in [Tf2] . In the present implementa- 
tion, the original proposal is extended to cover also initial-state splittings, finite parton masses, 
and QCD radiation off SUSY particles. 

Choices concerning the evolution parameter of the parton shower and the various scales 
entering running coupling constants, PDFs, etc. have been detailed, fixing the full algorithm. 
The kinematics of massive splittings has been presented in some detail, and the corresponding 
massless limits have been discussed. By direct comparison with some benchmark processes, 
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at first order in a s , the differences of the parton shower approximation with respect to exact 
results have been worked out. It has been shown that indeed the parton shower algorithm 
presented here reproduces the soft and collinear limits of the exact matrix elements and that 
differences between both results are non-singular terms only. Some first results with this new 
parton shower formulation have been presented and show very encouraging agreement with 
other models and with experimental data. 

In the near future, this new algorithm will be fully incorporated into the SHERPA frame- 
work and it will be made publicly available in the next releases of the code. This will also 
involve a more careful tuning of the shower parameters and the inputs of the hadronisation 
models provided by or linked to SHERPA, which surely will further improve the agreement with 
data. Planned is a detailed comparison against another new shower ansatz that is based on 
splitting colour dipoles [120] . and that is also being developed in the SHERPA framework at 
present. In addition, a full merging with multi-leg matrix elements in the spirit of [UJ will be 
implemented. It can furthermore be anticipated that this new shower implementation will lend 
itself to incorporation of MC@NLO-techniques 09]. 
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